IPliT-TlO/042 



Integrability and the AdS/CFT Correspondence 



Didina Serban 



InsUtut de Physique Theonque, CNRS- URA 2306 
C.E.A.-Saclay 
F-91191 Gif-sur-Yvette, France 

Abstract 



The description of gauge theories at strong couphng is one of the long-standing problems 
in theoretical physics. The idea of a relation between strongly coupled gauge theories 
and string theory was pioneered by 't Hooft, Wilson and Polyakov. A decade ago, Mal- 
dacena made this relation explicit by conjecturing the exact equivalence of a conformally 
invariant theory in four dimensions, the maximally supersymmetric Yang-Mills theory, 
with string theory in the AdS^ x background. Other examples of correspondence 
between a conformally invariant theory and string theory in an AdS background were 
discovered recently. The comparison of the two sides of the correspondence requires the 
use of non-perturbative methods. The discovery of integrable structures in gauge the- 
ory and string theory led to the conjecture that the two theories are integrable for any 
value of the coupling constant and that they share the same integrable structure defined 
non-perturbatively. The last eight years brought remarkable progress in identifying this 
solvable model and in explicitly solving the problem of computing the spectrum of confor- 
mal dimensions of the theory. The progress came from the identification of the dilatation 
operator with an integrable spin chain and from the study of the string sigma model. In 
this thesis, I present the evolution of the concept of integrability in the framework of the 
AdS/CFT correspondence and the the main results obtained using this approach. 



Contents 



1 Introduction: the gauge-string correspondence 4 

2 The A/" = 4 super Yang-Mills theory 11 

2.1 The action 11 

2.2 Symmetries 11 

2.3 Oscillator representation of PSU{2, 2|4) 13 

2.4 The dilatation operator as the Hamiltonian of an integrable spin chain . 14 

2.5 Bathe ansatz solution for the Heisenberg model 17 

2.6 The so(6) sector at one loop 20 

2.7 The full psu{2, 2|4) Bethe ansatz at one loop 21 

3 Perturbative integrability in = 4 SYM 24 

3.1 The notion of perturbative integrabihty and the asymptotic Bethe ansatz. 24 

3.2 The Inozemtsev model 26 

3.3 The three loop dilatation operator and Inozemtsev model 28 

4 Strings on AdSs x 29 

4.1 The Lagrangian 31 

4.2 Integrability of the classical sigma model 31 

5 Comparison between the gauge theory and the string results 34 

5.1 Strings in the plane- wave background. The BMN result 34 

5.2 Algebraic curves and the continuum limit of the Bethe ansatz equations: 

the su{2) principal chiral model 37 

6 All loop integrability in = 4 SYM 41 

6.1 The BDS conjecture 42 

6.2 BDS spin chain and the Hubbard model 45 

6.3 The dressing phase 49 

6.4 The all- loop asymptotic Bethe ansatz 50 

6.4.1 The conjecture 50 

6.4.2 The proof 52 

6.5 The dressing phase at all loops 56 

6.5.1 The quantization of the su{2) principal chiral model and the AFS 
dressing factor 57 

6.5.2 Higher loop corrections to the dressing phase 58 

6.5.3 The "magic" BES formula 61 

7 Strong coupling limit of the asymptotic Bethe An- satz equations 62 

7.1 The cusp anomalous dimension 62 

7.2 The generalized scaling function 65 

7.3 The strong coupling limit 66 

7.3.1 The resolvent 67 

7.3.2 The linearized BES equation 68 



2 



7.3.3 The functional equations 70 

7.3.4 The perturbative solution at strong couphng 72 

8 Finite size corrections and the thermodynamical Bethe ansatz 76 

8.1 Y system for the principal chiral field 78 

8.2 Y-system for the AdS/CFT Bethe ansatz 83 

8.2.1 Hirota equation and the Backlund transform for supersymmetric 

spin chains 84 

8.2.2 Hirota equation and the AdS/CFT integrable system 87 

8.2.3 Four loop Konishi operator and the Y-system 88 

9 Recent developments and open problems 90 
A Dynkin diagrams 91 
B Bethe ansatz equations as integral equations 92 



3 



1 Introduction: the gauge-string correspondence 



The relation between gauge theories and string theory is a long-standing and unsolved 
problem in theoretical physics. Understanding this relation will help us to handle gauge 
theories at strong coupling and consistently include gravity into the picture of quantum 
field theory. 

String theory itself was discovered in the attempt to describe the strong interactions, 
once the dual models [1] were proposed to incorporate the so-called Regge trajectories 
observed in experiments. The advent of the description of strong interaction in terms of 
Yang- Mills theory reduced the importance of string theory as a model of strong interac- 
tions, and for a while the interest in string theory stemmed mainly from the fact that it 
promised to give an unified description of all the forces in nature, including gravity. 

The idea of duality between strings and gauge theories originates in the work of 't 
Hooft [2], who understood that perturbative expansion of U{N) gauge theories in the 
large N limit can be reorganized in terms of a genus expansion of the surface spanned 
by the Feynman diagrams, expansion which is reminiscent of the perturbation expansion 
in string theory. The duality was also suggested by the strong coupling expansion of the 
lattice gauge theory of Wilson [3] and Polyakov [4] in terms of surfaces spanned by the 
color-electric fiuxes, which may end on quarks or make a propagating closed string. 

A simpler version of the gauge-string duality was extensively studied via the corre- 
spondence [5] [6] [7] between matrix models and non-critical strings, as well as through 
the study of the 2d quantum gravity. This subject is reviewed in [8] [9] [10]. Many of 
the ideas and techniques which were developed in this context are now transposed to the 
study of the modern, full-fiedged duality. 

The discovery of the D-branes as extended objects on which the strings can start 
or end [11] [12] gave a new impetus to the study of non-perturbative effects and string 
dualities and helped to understand the relation with black hole geometries and gauge 
theory. It was understood by Witten [13] that a stack of coincident D-branes is 
associated to a U{N) symmetric gauge theory. Shortly after, Klebanov [14] computed 
the absorption cross section of a stack of N D-branes, both from the point of view of a 
gauge theory and from the point of view of supergravity. For the case of D3-brane, the 
two computations gave exactly the same answer. 

Building on these ideas, Maldacena [15] conjectured in 1997 that the maximally 
supersymmetric gauge theory in four dimensions is exactly dual to type IIB strings in 
the AdS^ X background. He considered strings in a geometry created by a stack 
of coincident D3-branes, such that the theory on the branes is the four dimensional 
jV = 4 super Yang-Mills (SYM) theory. The geometry near the branes is that of an 
anti-de Sitter space times a five-sphere. The four-dimensional space can be recovered 
as the boundary of the AdS^ space, while the five-sphere is associated with the internal 
symmetry of the gauge fields. 

The conjecture by Maldacena confirmed the idea of Polyakov [16] [17] that the non- 
critical string theory describing gauge fields in four dimensions should be completed 
with an extra, Liouvillc-likc direction, which gives rise to a curved five-dimensional 
space. In order to insure the full reparametrization invariance of the Wilson loops, 
the four-dimensional gauge theory should be situated on the boundary of this curved 
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five-dimensional space. The idea that a five-dimensional theory is determined by a four- 
dimensional theory on the boundary rejoins an earlier proposal by 't Hooft, made in the 

context of the black hole physics and subsequently known under the name of holgraphy 
[18]. A similar property is associated with the topological field theories as the Chern- 
Simons theory, used to describe the quantum Hall effect. 

The gauge theory is characterized by two parameters: the coupling constant gyn and 
and the number of colors N. In the large N limit, or planar limit, they can be combined 
into a single parameter, the 't Hooft coupling constant A = Gym^ ■ These parameters 
are related to the string parameters by [15] 

47r^s = qIm = ^/N , i?V«" = A , (1.1) 

with Qs being the string coupling constant, a' the string tension and R the radius both 
of AdS^ and of the five- sphere. The identification of the parameters can be done by 
inspecting the supergravity solution [19] associated to N coincident D3-branes. 

The quantum M = 4 SYM theory is invariant under dilatations, so that the Poincare 
symmetry is extended to the conformal symmetry in four dimensions, with symmetry 
group 5*0(4,2) ~ SU{2,2). The symmetry group acting on the four copies of super- 
symmetry generators, the so-called R-symmetry, is given by SU{4:) ~ 5*0(6). Together 
with the supersymmetry generators, the symmetry group of the Af — A SYM theory 
is extended to the graded Lie group SU{2,2\A), the Af = A superconformal group in 
four dimensions. The bosonic subgroup 50(4,2) (S) 50(6) is the isometry group of the 
AdS^ X string background. Metsacv and Tseytlin showed [20] that the type IIB string 
theory on the AdS^ x background can be formulated as a non-linear sigma model on 
the coset superspace 5[/(2, 2|4)/50(4, 1) ® 50(5), thus proving that the two theories 
share the same global symmetry. 

A special role in the correspondence is played by the dilatation operator, which 
is associated to the Hamiltonian in the radial quantization. The eigenvalues of the 
dilatation operator are the quantum dimensions of the local operators in the conformally 
invariant theory. The gauge-invariant local operators in the gauge theory are nothing 
but the traces over the color indices of products of the fundamental fields. Under the 
correspondence, they arc associated to the states of the string, and their conformal 
dimensions are associated to the string energy. The dilatation operator corresponds to a 
non-compact generator of the global symmetry group, which can get quantum corrections 
and depends on the coupling constant. The other charges are discretized and correspond 
to angular momenta for the rotation of the string either on the AdS^ space or on 5^. 

The conjecture [15] is supposed, in its strongest form, to hold for any number of colors 
A^. In practice, it is easier to study the planar limit, when the number of colors goes 
to infinity, while the 't Hooft coupling constant A is kept fixed. In this case, according 
to eq. (1.1), the string coupling constant vanishes and we are dealing with the weak 
formulation of the AdS/CFT correspondence which associates the planar limit of the 
Af = A SYM theory to non-interacting strings in the AdS^ x 5^ background. In this 
thesis, we will concentrate exclusively on the planar limit. 

As it can be guessed from equation (1.1), the Maldacena duality is a strong/weak 
coupling duality, associating the perturbative regime A — >■ in the Af — A SYM theory 
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to the strongly coupled regime of the Metsaev-Tseytlin sigma model. This is the reason 
the duality is not easy to prove, since the regimes accessible by perturbative methods in 
the two theories do not match. However, due to supersymmetry, some of the states do 
not change with the coupling constant. They are the so-called BPS states and they were 
the main subject of the early checks of the correspondence. 

In 2002 the situation changed drastically when a new limit considered by Berenstein, 
Maldacena and Nastase (BMN) [21] allowed to explore a regime away from the BPS states. 
These authors have studied the excitations of around a point-like string configuration 
which moves with a large angular momentum J along a great circle on . The geometry 
seen by this fast moving string is the so-called plane wave geometry [22] , and string theory 
in this background is exactly solvable, the excitations being massive bosons/fermions. 
The solution depends on the effective parameter A' = A/ J^, the BMN couphng constant, 
which can take an arbitrary value. 

The same year, the first integrable structures where discovered, almost simultane- 
ously, in perturbative gauge and string theory. The first mention in the literature con- 
cerning integrabihty in the context of the planar J\f — A SYM theory goes back to 
Lipatov [23]. However, it is the work of Minahan and Zarembo [24], who mapped the 
one-loop dilatation operator to an integrable spin chain, which opened a whole new field 
of activity. Shortly after this discovery, Bena, Polchinski and Roiban [25] proved that 
the Metsaev-Tseylin sigma model is classically integrable. When it was understood that 
integrabihty extends to all operators at one loop [26] and that it persists at two and 
three loop order in gauge theory [27], it was tempting to conjecture that integrabihty is 
an all-loop feature. The integrable structures discovered in perturbation on the two sides 
of the correspondence would be then manifestations of the same object, which can be 
defined non-perturbatively. Given the power of the constraints usually associated with 
integrabihty, it is conceivable that integrabihty may allow to find a complete solution 
both to the problem of quantizing free strings in a particular curved background and 
that of solving a non-trivial quantum field theory in four dimensions. 

Other examples of the AdS/CFT correspondence were discovered recently, relating 
the J\f = 6 super-conformal Chern-Simons theory in 2-|-l dimensions and string theory in 
AdS^ X CP^ [28], as well as the AdS2,/CFT2 correspondence relating the type HB strings 
in AdSy, X S3 X M4 backgrounds, where M4 = T4 or S3 x Si, to some two-dimensional 
conformal field theory [29]. These theories are probably also integrable in the large N 
limit [30]. Other candidates to integrable correspondences were identified in [31]. 

The scope of this thesis is to present the concepts developed during the last eight 
years in connection with integrabihty in the context of the AdS/CFT correspondence. 
The relation between spin chains and sigma models is not new in physics. It was known 
for more than thirty years that the low energy physics of a Heisenberg antiferromagnet is 
captured by a two dimensional sigma model with a topological term [32]. However, the 
situation in the AdS/CFT correspondence is different, since there should be an exact, 
one to one mapping between the states of a spin chain and that of a sigma model. This is 
presumably made possible by the huge amount of symmetry, in particular the supersym- 
metry, involved in this particular theory. The integrable model behind the AdS/CFT 
correspondence includes as particular cases many known integrable models, but it have 
features that go beyond any particular model which was already studied. It combines 
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features of super-symmetric, non-compact spin chains with those of non-relativistic inte- 
grable sigma models, and it may serve to generahze and unify these different approaches 
to integrable models. The development of the subject itself came from the interplay of 
several approaches and benefitted from experts coming from various backgrounds and 
with different points of view. Crucial information was obtained at weak coupling in 
the gauge theory and in string theory from state-of-the-art perturbative computations. 
Notions and techniques from integrable spin chains, including super-spin chains, two di- 
mensional sigma models and matrix models were combined together and helped to build 
a coherent picture of the model. 

Given the fact that there is no information available for finite coupling constant 
neither from the gauge theory nor from string theory, the form of the equations describing 
the spectrum of conformal dimension was obtained, at least in the beginning, by trial 
and errors and by conjectures and checks rather than by constructive proofs. The formal 
proof of integrability is still lacking, but the fact that the same set of equation is able 
to reproduce both weak and strong coupling results is remarkable and cannot be an 
accident. 

After the conjecture of Beisert, Kristjansen and Staudacher [27] that the dilatation 
operator of = 4 SYM is integrable at all loop, the search of the all-loop equations de- 
termining the spectrum has started. The dilatation operator at higher loops corresponds 
to a long range Hamiltonian, with the range of the interaction growing with the loop 
order. The closest known example of this type is the Inozemtsev Hamiltonian [33]. This 
system can be diagonalized using Bethe ansatz-like equations only for infinite length, 
with corrections which are exponential in the system size. These equations will be called 
in the following asymptotic Bethe ansatz equations. Serban and Staudacher showed [34] 
that the dilatation operator can be approximated, up to three loops and for a partic- 
ular class of operators, by a combination of the integrals of motion of the Inozemtsev 
model and therefore could be diagonalized by the Inozemtsev asymptotic Bethe ansatz. 
The first corrections in the BMN coupling constant A' coincided with those obtained 
from string theory [35, 36, 37], while the third loop order disagree with the string result 
due to order-of-limit problem. Staudacher [38] emphasized the link between the asymp- 
totic Bethe ansatz and the concept of scattering matrix for magnons, concept which 
subsequently proved to be very fruitful. In particular, using the S matrix allowed to 
circumvent the difficulties related to the changing length [39] which renders the usual 
algebraic Bethe ansatz approach to spin chains inoperable in the present context. 

At the same time, Kazakov, Marshakov, Minahan and Zarembo [40] studied the so- 
called finite gap classical solutions for the su{2) principal chiral model and expressed 
them in the elegant language of algebraic curves. Since the continuum limit of the Bethe 
ansatz equations leads also to algebraic curves, it became possible to compare directly 
the sigma model predictions with the spin chain predictions. The same method was used 
for the s/(2) sector [41] and later for the whole psu{2,2\A) algebra [42]. 

The finite size corrections in the language of the spin chains correspond to loop 
corrections for the string sigma model. These were extensively studied [43, 44, 45, 46, 
47, 48, 49], and shown to coincide, at least for the first few orders in the BMN coupling 
constant, where the order-of-limit problem does no show up yet. 

The direct comparison between the algebraic curves of the sigma model and the 
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Bethe equations made possible to conjecture the all-loop Bethe equations. First, Beisert, 
Dippel and Staudacher (BDS) [50] were able to conjecture the dispersion relation and a 
set of Bethe ansatz equations which reproduce succesfully some of the all-loop features. 
However, it became soon clear [51] that a part was missing from the BDS conjecture, the 
so-called dressing factor. The BDS equations were soon extended to all sectors by Beisert 
and Staudacher [52], and the Beisert-Staudacher equations were rederived by Beisert 
[53] by the assumption that the excitations of the spin chain are subject to a centrally 
extended su(2\2) x su(2\2) symmetry. Subsequently it was found [54] [55] that these 
equations bear deep links, not yet entirely clarified, with the one-dimension Hubbard 
model. The BDS dispersion relation was confirmed, from the string side, by Hofman 
and Maldacena [56], (see also [57]), who coined the term giant magnon for the magnon 
excitations with momentum of the order of unity. The giant magnons correspond to 
string excitations were the effect of the curvature of the sphere is manifest. Arutyunov, 
Frolov and Zamaklar [58] showed that the Bethe ansatz equations for the strings should 
be essentially the same as that for the spin chain, since the symmetry of the excitations 
in the string sigma model is the same as that of the magnons in the spin chain [59] . 

The dressing phase was first determined at leading order at strong coupling by Aru- 
tyunov, Frolov and Staudacher [51]. A partial answer for the one-loop string corrections 
were obtained by Beisert and Tseytlin [48], building on previous work [44], and the full 
one-loop answer was obtained by Hernandez and Lopez [60]. Janik [61] used the su{2\l) 
symmetry and the quantum group structure to derive the equivalent of the crossing 
equation obeyed by the dressing factor. A solution to the crossing equation was found 
by Beisert, Hernandez and Lopez (BHL) [62]; this solution was shown by Beisert, Eden 
and Staudacher (BES) [63] to obey the right structure at weak coupling. The BHL/BES 
solution of the crossing equation allowed to check the fourth loop order in the cusp 
anomalous dimension, obtained by Bern et al. [64] and passed a number of other consis- 
tency checks. Dorey, Hofman and Maldacena [65] checked that the BHL/BES solution 
has the expected set of singularities and provided an useful integral representation. Very 
recently, Arutyunov and Frolov [66] verified that the BHL/BES solution satisfies the 
crossing equation and Volin [67] showed, conversely, that the solution of the crossing 
equation with the minimal number of singularities in the physical strip is the BHL/BES 
solution. Finding the minimal solution of the crossing equation provided the interpo- 
lation between the weak coupling and strong coupling, at least for spin chains of large 
length. 

A considerable effort was spend to analyze the Bethe ansatz equations at strong 
coupling, especially concerning the cusp anomalous dimension^. An integral equation, 
the BES equation, giving the leading of the cusp anomalous dimension part at large 
magnon number (or Lorentz spin) was written down by Eden and Staudacher [70] and 
corrected with the dressing factor by Beisert, Eden and Staudacher [63]. This equation 
was analyzed numerically by Benna et al. [71], who confirmed the first two orders 
obtained previously with string theory techniques [72] [73] and gave a prediction for the 

^The cusp anomalous dimension governs the renormalization of Wilson loops with cusps [68] and 
controls the infrared behavior of the scattering amplitudes in gauge theories [69] . They are ubiquitous 
observables in gauge theories, appearing among others in the logarithmic scaling of the anomalous 
dimensions of operators with high Lorenz spin. 
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third order which was found later [74, 75]. Analytical results were much more difficult 
to obtain. The first order was reproduced by a series of works [76, 77, 78, 79], the next 
term was obtained by Casteill and Kristjansen [80] and Belitsky [81], and the third term 
and a recursive procedure for the next terms by Basso, Korchemsky and Kotanski [82] 
and in [83]. Another object extensively studied at strong coupling, from the point of 
view of the Bethe ansatz and the string theory, is the generalized scaling function which 
also appears in the study of the logarithmic scahng in the sl{2) sector of the theory. An 
equation giving the generalized scaling function was written down by Preyhult, Rej and 
Staudacher (FRS) [84], and rederivcd and the series expansion in the extra parameter 
j was studied in [85]. Basso and Korchemsky [86] showed that at small values of the 
parameter j the FRS equation reduces to the equation for the energy of the 0(6) sigma 
model, confirming the relation with the 0{6) model identified previously by Alday and 
Maldacena [87]. At large values of the parameter j, of the order of the coupling constant 
g, the behavior of the generalized scaling function was computed from string theory by 
Roiban and Tseytlin [88]. The same quantity was computed form the Bethe ansatz, 
using different approaches, by Gromov [89], Bajnok et al. [90] and Volin [91], who found 
a mismatch in one of the coefficients with the string computation [88], revisited and 
corrected in [92]. 

Among the most recent developments are the study of finite size corrections and 
the relation between gluon amplitudes and integrability. The finite size corrections in 
AdS/CFT were considered first by Ambjorn, Janik and Kristjansen, [93], from the point 
of view of an integrable field theory. Janik and Lukowsky [94] proposed to use the method 
initiated by Liischcr [95, 96] to compute the finite size corrections. This approach was 
successfully carried out for Konishi operator at four loops [97] and at five loops [98]. The 
four loop result coincides coincides with the direct diagrammatic computations [99, 100], 
while there is no five loop direct computation in the gauge theory. 

Another approach for computing finite size corrections is connected with the Thermo- 
dynamic Bethe Ansatz, and the underlying idea put forward by Alexei Zamolodchikov 
[101] is to exchange the role of the finite temperature with that of a finite length by 
a transformation similar to the modular transformation in 2d CFTs. In the context of 
AdS/CFT, this approach was suggested already in [93] and a systematic implementation 
was initiated by Arutyunov and Frolov [102]. Gromov, Kazakov and Vieira [103] showed 
for the case of the su{2) principal chiral model that these equations, named the Y-system, 
can be used to determine the finite size corrections. When applied to AdS/CFT, this 
method allowed to reproduce the four loop [104] and five loop [105] [106] corrections for 
the Konishi operator previously obtained by the Liischer method. 

While most of the present developments related to integrability concern the deter- 
mination of the conformal dimensions and the study of the dilatation operator as an 
integrable spin chain, there are signs that integrability has deeper consequences on the 
structure of the supersymmetric J\f = i gauge theory. It was discovered [107] that the 
four-gluon amplitude has a particular structure, in which the one-loop contribution ex- 
ponentiates [108]. This property, checked at strong coupling by Alday and Maldacena 
[109], docs not hold for six or more gluons [110, 111]. Sokatchev and his collaborators 
[112] discovered that the integrals which enter the expression of the amplitudes possess 
a dual conformal invariance. Shorthly after, it was discovered that the existence of this 
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dual conformal invariance is the consequence of a duality between the amplitudes and 
the Wilson loops constructed from light-hke segments with cusps. This explains the 
appearance of the cusp anomalous dimension in the expression of the amplitudes given 
in [108]. At strong coupling, the relation between amplitudes and the Wilson loops with 
cusps appear via a T-duality [109], which was extended to a fermionic T-duality in [113]. 
The action of this duahty on the integrable structure of the AdS^ x superstring was 
analyzed in [114] and reviewed in [115]. 

Drummond, Henn and Plefka [116] proposed that the generators of the super-confor- 
mal symmetry, together with those of the dual super- conformal theory, constitute the 
generators of a Yangian. This would imply the existence of an integrable structure of 
the amplitudes. An integrable model for the high energy limit of the amplitudes was 
proposed by Lipatov [117]. 

Non-planar corrections to the dilatation operator, corresponding to interactions split- 
ting and joining the spin chains, were investigated in [118] [119]. 

The material is structured as follows: in section 2, the basic facts about the J\f = 4 
super Yang-Mills theory are presented, with an emphasis on the symmetry and a simple 
realization of the symmetry in terms of oscillators. Section 2 contains also a presentation 
of the dilatation operator at one loop, together with the exact solution for its spectrum 
in terms as the spectrum of an integrable spin chain. Section 3 contains the definition 
of perturbative integrability for the dilatation operator and the comparison with the 
Inozemtsev model. Section 4 is devoted to strings on the AdS^ x S^, its definition as a 
coset sigma model and the classical integrable structure, including the finite gap solution 
in terms of algebraic curves. In section 5 are presented the early results concerning the 
comparison between spin chain solutions and the string solutions for strings on plane- 
wave background, as well as the comparison between the Bethe ansatz equations for 
spin chains in the su{2) sector and the finite gap solutions for the su{2) principal chiral 
model. In section 6, the main results concerning the asymptotic all-loop Bethe equations 
are described, together with the solution for the dressing phase. Section 7 is devoted 
to the analysis of the Bethe ansatz equation at strong coupling, concerning the cusp 
anomalous dimension and the generalized scaling function. In section 8 we present the 
TBA approach for the su{2) principal chiral model and then for AdS/CFT, including the 
derivation of the wrapping correction for the Konishi operator at four loops. In section 9 
we briefly discuss the recent developments concerning the finite-size corrections and the 
remarkable properties of the amplitudes. 

In order to keep the text short we have skipped some of the details. For a more in- 
depth presentation of different aspects, one can consult the original works or the available 
reviews, e.g. the reviews by Plefka [120] [121], the PhD thesis of Beisert [122] or the 
recent review on superstrings in AdS^ x by Arutyunov and Frolov [123], as well as 
the special volume dedicated to the subject by the journal J. Phys A [124]. 
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2 The J\f = 4 super Yang-Mills theory 

2.1 The action 

The action for the Af — A gauge theory in four dimensions can be obtained by dimensional 
reduction from the ten dimensional Af — 1 gauge theory 

S^l AQTrFMivF^^ + ^IVVr^^^MV') (2.1) 

where the covariant derivative is defined as 

T^M^dM- igyu [Am, ] (2.2) 

and is a sixteen-component Majorana-Weyl spinor in ten dimensions. Upon dimen- 
sional reduction, six of ten components of the gauge field become scalars, while the 
sixteen-dimensional spinor decomposes into four copies of left and right two-component 
spinors in four dimensions 

Am, M = 1,...,4 ^ A^, /i = 0,...,3 (2.3) 
Am, M = 5,...,10 ^ i = l,...,6 

ipA , A=l,...,16 ij'^, iJa,a , a=l,...,4, a,a = 1,2 

The action for the four dimensional gauge theory becomes 

d^x (-Tr Fi.^F'"' + -Tr - ^Tr $j] [$\ (2.4) 

\4 2 4 

+Tr^V^I)^^, - !^Tr<Vj$\^,] - !^Tr a^.^'^f^,, ^'']) 

where a'^ and a* are the chiral projections of the gamma matrices in four and six dimen- 
sions respectively^. The fields A^, $i and ipa form a supermultiplet. 

2.2 Symmetries 

The action (2.4) is Poincare invariant and scale invariant. The classical dimensions of 
the fields are 

K] = [*,] = i m^l- (2.5) 

The theory remains scale invariant upon quantization, the beta function being zero, 
and the theory is conformally invariant at quantum level. The conformal group in four 
dimensions is 5*0(4,2) ~ SU{2,2). This group contains two SU{2) components, with 
generators and L^, the dilatation operator, D, the translations, and the special 
conformal transformations, K^. The internal symmetry which rotates the six scalars into 
one another, or the R symmetry, is given by SO{6) ~ 5'C/(4). Taking into account the 

^For a list of properties of the chiral projections of the gamma matrices, see [125]. 
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super-translations Q^, Qaa and the super- version of the special conformal transforma- 
tions, S^, we obtain the total symmetry group as P5'C/(2, 2|4). The structure of 
the generators can be schematically represented as 

L, I, P, K, D Q, S\ 

Q, S R J ^ ' 

where the generators in the diagonal blocks are bosonic and the ones in the anti-diagonal 
blocks are fermionic. The generators have a definite dimension, which is not modified 
by the radiative corrections 

[D] = [L] = [L] = [i?] = , [P] = 1 , [K] = -1 , [Q] = 1/2 , [S] = -1/2 . (2.7) 

It will be useful in the following to represent the the fields in the spinorial notation, 
using the properties of the the chiral projections of the gamma matrices 

V^ = a^^V^f,, $, = af$„b, etc (2.8) 

Conformal dimensions and the dilatation operator 

We are going to consider gauge invariant, local operators, which are traces over the 
gauge group of products of operators, e.g. 

Oi,^i,...aM = Tr[^,,{x)V^^,,{x)...^/Ja{x)^^Ax)] , (2.9) 

Multiple trace operators occur as well, but in the planar limit we can restrain ourselves 
to the single trace operators. 

The operators organize in super-multiplets of the PSU{2,2\4:) symmetry. The op- 
erator with the lowest dimension in the multiplet is called a superconformal primary 
operator. 

The unitary representations of the superconformal algebra are labeled by the quantum 
numbers of the bosonic subgroup 

SO{3, 1) X SO{l, 1) X SU{A) (2.10) 
(si,S2) A(^) [ri,r2,r3] 

with si and S2 being half-integer, A{g) is the positive conformal dimension and [ri, r2, r^] 
are the (integer) Dynkin labels of SU{A). 

The unitary representations of the superconformal group were classified by Dobrev 
and Petkova [126] [127]. There are three discrete series of representations, for which at 
least one of the generators Q commutes with the superconformal primary operator. These 
are the BPS or short multiplets, and their dimension does not vary with the coupling 
constant. They correspond to protected operators. A fourth series of representations is 
continuous, and their dimension vary with the coupling constant. These are the non- 
BPS operators, and computing A{g) is the main concern of the activity reported here. 
A(0) is the classical dimension of the field, which is the sum of classical dimensions of 
the components of a composite field. 
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2.3 Oscillator representation of PSU{2, 2|4) 

It is sometimes useful to employ an oscillator realization of PSU{2,2\A) [128], which is 
valid for the free theory, qym = 0. Consider two sets of bosonic oscillators, (a", aj^) and 
{W, ftt ), with a,a— 1,2 and a set of fermionic oscillators, (c", c^) with a = 1, 2, 3, 4, with 
standard commutation relations 

[a", a\] = 5-^ , r , 6t ] = S'^ , {c«, cj} = 5^ . (2.11) 

The fundamental representation of S'f/(2,2|4) and its conjugate can be realized in the 
space spanned by the eight oscillators {A"^) = {a,b^,c) and (^^4) = {a'^,—b,c^) respec- 
tively. The particle-hole transformation on the b oscillators is related to the non-compact 
nature of the group. The generators of the algebra are, in this representation, given by 
the bihnears A^^A^ and they act on the states by commutator. The su{2) x su{2) x su{A) 
part is realized simply as 

Li = aia^-^ayS^^, (2.12) 
Li = bib^-lb\b^6i, (2.13) 
K = cy-\cy5l. (2.14) 

The dilatation operator is given by 

Do^l + ^{na + n,) , (2.15) 

while the central charge, which should be zero for PSU{2, 2|4), is given by 

C = l-^(n„-n6 + ne) , (2.16) 

which means that for the representation we consider, we have Uc -\- Ua — ni, — 2 . The 
off-diagonal generators are 

Pap = <^ib\. K<^^ = a"b^, (2.17) 
Ql - alc\ S^ = cla\ (2.18) 
Q^a = b\ci, S'^^^b'^c^. (2.19) 

On this representation, it is easy to check the commutation relations of the psu{2, 2|4) 
algebra 

[S:, PJ = Q^„, [K-^ Q^a] = 5? S: (2.20) 

pj = 5| [K^^, g^] = (2.21) 

{5f , S^'^} = 5^, K^^ {Q' Q^a} = 51 Pp^ (2.22) 
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and 

\K-^, Pp^] = 5yLi + 5i + 5"^5iD (2.23) 

{St. Q'p) = + 5tR'a+ D (2.24) 

J = 51 Li -5iRl + Ut 6i D (2.25) 

The different A/" = 4 operators can be realized as states in an infinite dimensional 
PSU{2,2\4) representation as follows: 

a) the bosons Z = $i + i$2, ^ = ^3 + and y = $5 + i^g by states with two 
fermions 

z~44|o), x~44|o), y~44|o). 

b) the fermions 

c) the covariant derivatives 

One can make a particle-hole transformation on the fermions and c"^ 

4 ^ c^ 4 = c^ (2.26) 

such that the highest weight state, associated to the operator Z, is the "empty" state, 
annihilated by all the annihilation operators a*, 6*, c*, rf* with i = 1,2. The reference 
state Z is manifestly invariant under su{2) x sm(2), and in fact it is invariant under the 
action of s-u(2|2) x sm(2|2) generated by the oscillators (a, c) and {h,d) respectively. 

Dynkin diagrams and the oscillator representation 

There is a close relation between the sl{N\M) root system, as described in Appendix 
A and the oscillator representation. 

For each e {6) direction defined in Appendix A one can choose a bosonic (fermionic) os- 
cillator. The bosonic roots will be then represented by bilinears in bosons or of fermions, 
while the fermionic root will contain one boson and one fermion. Depending on the real 
form one wants to choose for the bosonic part, one can make a particle- hole transfor- 
mation for a set of the bosons. For sm(2,2|4) one recovers the construction from the 
previous section. 

Among the different Dynkin diagrams for sii(2,2|4), the two shown in the figure 1 
will be especially useful for writing the Bethe ansatz equation. They show the explicit 
decomposition with respect to sm(2|2) sm(2|2), algebra which can be obtained upon 
deleting the middle node of the two diagrams. 

2.4 The dilatation operator as the Hamiltonian of an integrable 
spin chain 

The basis of local operators can be written as a collection of words of arbitrary length 
made by the letters corresponding to basic operators 

Tr {ZZXV3Yijjiij:2 . . . VqXY) {x) . (2.27) 
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Figure 1: Two Dynkin diagrams for psu{2, 2|4). 

taken at a particular point x in the four-dimensional space-time. Each basic operator can 
be interpreted as spin state associated to a site of a periodic spin chain. The dilatation 
operator D maps these spin states into one another. The eigenvalues of the dilatation 
operator are the conformal dimensions A^{g) of the operators O'^lx), i.e. 

mx)0\y)) ~ ^3^^ (2.28) 

The appearance of non-trivial dimensions is a consequence of renormalization. Although 
the beta function of the theory is zero, there are singularities at higher loops, which need 
to be removed by renormalization. For example, the correlation function (2.28) can be 
computed in perturbation theory as 

{0\x)0%)) ~ (1 - 2/A^ln(A|x - y|) + . . .) (2.29) 

\X — y| 

such that the perturbative expansion of the conformal dimension is given by 

^\g) = + g^^^ -r... = ^^ + 5^^{g). (2.30) 

Here we had to introduce the scale A which plays the role of the UV cutoff. Alterna- 
tively, the anomalous part of dilatation operator SD{g) can be read off the wave function 
renormalization for the operators 

Of^,{x) = Z% C)^(x) = [e&D{9)^-KYj, C)^(x) (2.31) 

The dilatation operator can be therefore computed in perturbation and it can be 
interpreted as a spin Hamiltonian which maps a spin state into another spin state. 
Because of the cyclicity of the trace, only the states which are invariant by translation 
should be considered. The problem of computing the (anomalous) dimensions is then 
reduced to diagonalizing the Hamiltonian of the spin chain. 

Let us remind that the dimensions of the symmetry currents J{(j) should not be 
modified in perturbation theory, so that 

\D{g), J{g)\ = [Do, J{g)\ = dim(J) J(^) , (2.32) 

where is the classical part of the dilatation operator and the dimensions dim(J) 
are given by equation (2.7). This means in particular that the anomalous part of the 
dilatation operator -D(g') — -Do commutes with all the generators of psu{2, 2|4) 

[D{g)-Do,J{g)]=0 (2.33) 
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and therefore it can be considered as a psu{2, 2|4) symmetric Hamiltonian. In the same 
time, D(g) is an element of the psii(2,2|4) algebra, a feature that distinguishes it from 
the usual spin chains. 

Another consequence of equation (2.33) is that the one-loop part of the dilatation 
operator is invariant under the classical part of the superconformal symmetry 

p2, Jo] = . (2.34) 

We notice that a change of sign g — ^ —g can be absorbed in the A/" = 4 action by 
a sign change of the scalar fields $j — ?> — Such a change should not modify the 
conformal dimensions ^{g) , so we conclude that they depend only on g'^. 

Hg)^Y.9"'^2k. (2.35) 

fc>0 

On the other hand, the dilatation operator itself may contain odd powers of g. 

The mixing matrix at one loop was first computed in [21], and it was noticed by 
Minahan and Zarembo [24] that, for the so(6) sector spanned by the 6 scalar bosons, 
the corresponding Hamiltonian is the integrablc so(6) spin chain, solved by Reshetikhin 
[129] [130]. The full one-loop Hamiltonian was found by Beisert [26] and the full one-loop 
Bethe ansatz was written in [131]. 

The mixing matrix and the sectors 

The spin chain can be represented in terms of oscillators, with a copy of oscillators 
at each site of the spin chain. While at = there is no interaction between sites and 
the dilatation operator is just the sum classical dimensions at each site (2.15), when the 
interaction is turned on, the dilatation operator has a non-trivial expression in terms 
of the creation/annihilation operators. Some restrictions may be imposed on the total 
oscillator content of the states, and some of these restrictions are conserved by the full 
dilatation operator. The dilatation operator is a Cartan generator for the ps-u(2,2|4) 
algebra, and therefore its action cannot change the ]9sm(2,2|4) charges of the states 
on which it acts, charges that are related to the oscillator content of the states. The 
restrictions give rise to sectors which are closed in perturbation theory. Let us remind 
that any state in the spin chain has to satisfy the constraint 

ria -\- Uc — ni, — 2L . (2.36) 

• The BPS states: in the particular case 

= nci = L , =nh^=0 (2.37) 

the constraint (2.36) is saturated and the action of the anomalous part of the 
dilatation operator is zero. They correspond to states \ZZZZ...ZZZ). These are 
the BPS states protected by supersymmetry, together with any state obtained from 
it by a,SO{6) rotation. 

• The su{2) sector: 

na, =nh^=0 and Uc,^ = n^^ + n^^ = L (2.38) 
These states are of the type \ZXXZ...ZZZ) . 
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• The sl{2) sector: 

n-ai = ribi = M, = rib^ = 0, and rics = = L; , (2.39) 
with states are of the type \V^ZZ...ZZZ). 

• The sector: 

= ribi = and = ria^ + = L. (2.40) 
The states are of the type \ZiIjiiIjiZ...Z Z Z) with ■^i = ^31. 

• The sti(2|l) sector: 

n-ai = rifti = and ric^ = riaj + ric^ +nc2^ L (2-41) 
spanned by states of the type \ZipiipiZ...ZXZ). 

• The su{3\2) sector: 

= and ricg = + + n^^ + n^.^ + n^^ ^ L (2.42) 

spanned by states of the type \Zil)i'il)2Z...ZXY Z) with ipi = ip^i and 1^2 = ip32- The 
transformation \...XYZ...) \...'ijji'ijj2---) aUowed by diagrammatics corresponds to 
2- This transformation conserves the psu{2, 2|4) charges but reduces 
the length by one unit, L — > L — 1. This sector is the smallest sector with variable 
length. 

• The psm(1, 1|2) sector: 

^nb^^O and ricg - Ub^ + ric^ + = L . (2.43) 

The states are of the type |X>^ZV'i'0i^...^XZ) with t/Si = -i/Sj. 



2.5 Bethe ansatz solution for the Heisenberg model 

When reduced to the operators containing only the two complex combinations of the 
scalar bosons, for example Z and X, the one-loop dilatation operator reduces to the 
Heisenberg Hamiltonian 

//r^^) = 2/X;(l-n,+0, (2.44) 
1=1 

where 1 is the identity operator and Pi^i+i permutes the states at sites / and / + 1. 

The solution for the Hamiltonian (2.44) can be found by the so-called coordinate 
Bethe ansatz, which was originally employed by Bethe [132]. One starts with a reference 
state, where all the spins are in the "up" state | t). In terms of the J\f — A SYM, one 
starts with a trace composed uniquely from Z operators. 

\n) = I nn • • • t) (2.45) 
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Starting with this state, one can generate the whole Hilbert space by reversing an arbi- 
trary number of spins. A reversed spin 



i;) = ^"it) 

will be called a magnon. Since the Hamiltonian (2.44) conserves the number of reversed 
spins, it can be diagonalized on a space with fixed number of magnons. Consider first 
the one-magnon states 

L 

*(p) = ^e'^V,-|Q) (2.46) 

k=l 

which are eigenstates of (2.44) if p = 27rn/L, n e Z, with eigenvalue 

E(p) = 2/(2 - e'P - e-'P) = 8/ . (2.47) 

For a state with two magnons, the eigenfunctions are found with the following ansatz 
*(Pi,P2) = Yl [^(Pi,P2)e'^^i*^i+^^*^^) +>l(p2,Pi)e'^*'^'=i+^i*^^)] a;^a^jn){2AS) 

l<ki<k2<L 
l<k2<ki<L 

This ansatz is clearly valid when ki <^ k2, since the action of the hamiltonian (2.44) 

couples only two nearby sites. The condition that it is valid everywhere determines the 
ratio of the two coefficients A{p2iP\) and A{pi,p2) which defines the scattering matrix 
of two magnons 

, _ A(p2,Pi) e^P^+^^^ - 2^P- + 1 
^^^^'^^^ = A^) = + 1 ^2.49) 

while the invariance under translations which is implicit on the summation imposes that 

e'v^L ^ s{p2,pi) , e'f^^ = S{pi,p2) . (2.50) 

The equations (2.50) are valid for any spin model with local interaction. What is remar- 
cable about the model (2.44) is that an ansatz of the type (2.48) is valid for any number 
M of magnons. The wave function in a different chamber /i;p(i) < /cp(2) < . . . < kpi^M) 
can be obtained from the one in /ci < /c2 < . . . < ku just by the inverse permutation 
of the momenta Pi,P2, ■ ■ -Pm Pp-^{i),Pp-^{2), ■ ■ ■Pp-'^{m)- Since any permutation P 
can be written as a product of transpositions, the elementary object is the two magnon 
scattering matrix (2.49). The magnon momenta are determined by the equations 

M 

e'Pr.L^ -Q S{pn,Pm), n = l,...M (2.51) 

The property of factorizcd scattering is related, at least for the systems with a finite 
number of degrees of freedom, to the existence of a number of conserved quantities equal 
to the number of degrees of freedom, that is, to integrability. 
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A different way to obtain the Bethe ansatz equations (2.51), as well as the conserved 
quantities, is via the algebraic Bethe ansatz. For a pedagogical review of different aspects 
of the algebraic Bethe ansatz, sec [133]. One starts with a matrix R(u) acting onV<S>V, 
with V = <C? which satisfies the Yang-Baxter equation 

Rr2{u - v)R,3(u)R23(v) = R23(v)Ri3(u)Ri2(u - v) (2.52) 

where the subscript indicates the space where the matrix act. In the su{2) case, such a 
solution is provided by 

Ri2{u)^u + iPi2 (2.53) 

where P12 is the permutation operator. For u — ±i the matrix R{u) is proportional to a 
projector. Starting with the R matrix, one can construct the monodromy matrix To(it) 

by 

To{u) = U{u)L2{u) . . . Ll{u) , with L„(m) = Ron{u) (2.54) 
which satisfies the Yang-Baxter equation in the form 

Roo'iu - v) To{u) To/(f ) = To>{v) To(m) Roo'iu - v). (2.55) 

A consequence of the equation (2.55) is that the trace of the monodromy matrix, called 
the transfer matrix, commute with itself for different values of the spectral parameter u 

[T{u),T{v)] = , T{u) = TroToiu) , (2.56) 

such that the transfer matrix T{u) generates the conserved charges. A special point is 
u — 0, where the the conserved charges generated by the logarithm of transfer matrix 
are local, e.g. 

^(L - H2) = InT(^/)U=o . (2.57) 

The common eigenvectors of the conserved quantities can be constructed from the ref- 
erence vector 1^2) as follows. The monodromy matrix can be written as 

^»(-) - ( ) 

The conserved quantities are generated by 

T{u) = A{u) + D{u) (2.59) 

while the eigenvectors are build as 

M2, . . . um) = B{ui)B{u2) . . . BiuMm . (2.60) 

The vectors (2.60) are eigenvectors of the operators (2.59) if the rapidities Ui, U2, ■ ■ ■ um 
satisfy the equations 

Uk + 1/2 \ _ -|-|- Uk-ui + i (261) 
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which are identical with the set of equations (2.51) upon the identification 

= ^ cot I . (2.62) 
The energy is given, in these variables, by 

M 

£;f(2)^2/y ^-i— . (2.63) 
«fe + 1/4 

2.6 The so(6) sector at one loop 

Minahan and Zarembo computed the one loop dilatation operator for all traces involving 
operators in the so(6) sector^, that is, involving any of the bosonic operators $j with 
i = 1,...,6, 

!¥($,,$,,...$, J (x). (2.64) 

The dilatation operator acts as a Hamiltonian for a so(6) spin chain with the spins in 
the fundamental representation, 

L 

^soi6) ^g2j2 (Ki,i^, + 2 - 2P,,,+0 (2.65) 
1=1 

where the trace operator Ki^i^i contracts the indices on the nearby sites I and I + 1 
while the permutation operator permutes the indices at the sites / and / + 1 

Pit:: - ^t'^iiu ■ (2-67) 

This spin chain is integrable and it solution was obtained by Reshetikhin [130]. It also 
can be obtained [134] by the method of nested Bethe ansatz, proposed by Kulish and 
Reshetikhin [135]. The algebraic Bethe ansatz method from the previous section can be 
extended to algebras of higher rank and it associates a set of creation operators B^'^\ur) 
with each node r of the Dynkin diagram. It is called nested because at each step of the 
procedure the rank of the algebra is reduced. The Hamiltonian (2.65) is diagonalized in 
terms of the rapidities satisfying the following system of nested Bethe ansatz equations 

Ml , . M2 

U2,k + l/2 \ ^ TT- U2,k - U2,l + I Y\ U2,k - Ul,l - V2 -pi- 

U2,k-i/2j ~ ' -„,.^./9 11 





- Ul,l 


+ 1 


Ul,k 


- Ui^i 


— i 


U2,k 


- U2,l 


+ i 


U2,k 


- U2,l 


— i 


U3,k 


- U3,l 


+ i 


U3,k 


- U3,l 


— i 



M3 , . M2 





- U2,i - i/2 


Ul,k 


- U2,i + i/2 


U2,k 


- ui^i - i/2 


U2,k 


- ui^i + i/2 


U3,k 


- U2,i - i/2 


U3,k 


- U2,i + i/2 



^At higher loop order this sector is not closed and it mixes with sectors containing fermions. 
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1 

Figure 2: The Dynkin diagram and the Dynkin labels ct^ • w corresponding to the so(6) 
Bethe ansatz (2.68). 

The eigenvalues for the energy are given by 

Ma ^ 
i=l ' 

The su{2) Bethe equations can be obtained from (2.68) by restricting Mi = M3 = and 
M2 — M. Due to the invariance of the trace which defines the local operators in the gauge 
theory, one has to retain only the states of the spin chain with zero total momentum, 
which are therefore invariant by translation. This condition can be reformulated as 

Ma , • /o 

and has to be supplemented to the equations (2.68). 

The structure of the nested Bethe ansatz (2.68) is closely related to the data from 
the Dynkin diagram. The equations were extended to any simple Lie algebra in any 
representation by Ogievetsky and Wiegmann [136] . For any simple root in the Dynkin 
diagram there is a set of rapidities Ur^ki with k — 1 . . . M^. These rapidities obey the set 
of equations 

Ur,k + io^r ■ w/2 \^ ~ TT rf ^'''^ ~ " """'^^ (2 71) 

where w is the vector characterizing the representation, encoding the Dynkin labels of 
the highest weight vector. The corresponding energy is given, up to a multiplicative and 
additive constant, by 

Mr / . .X 

E-yy — - — 1 — rV (2-72) 



r k=l 



The Bethe equation (2.68) can be alternatively seen as equations for vector repre- 
sentation so(6) with simple roots €1 — 62, ^2 — ^3, ^2 + ^3 and weight w — ei, or for the 

antisymmetric representation of su{4) with roots ei — 62, 62 — 63, 63 — 64 and weight 
w = ei + 62- The corresponding Dynkin diagram is represented in figure 2. 

2.7 The full psu(2,2\4:) Bethe ansatz at one loop 

The action of the one loop dilatation operator for the whole psu{2, 2|4) states was com- 
puted by Beisert [26]. The main ingredient is the computation of the one- loop dilatation 
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operator in the sl{2) sector, which then can be hfted to an action on the whole psu{2, 2|4). 
In terms of the J\f — A fields, the sl{2) sector corresponds to one of the complex scalar 
fields, say Z, acted upon by an arbitrary numbers of covariant derivatives V. The corre- 
sponding representation has weight —1 or spin s = —1/2, and it is a infinite dimensional 
heighest weight representation of sl{2). The sl{2) Hamiltonian was determined by Beis- 
ert, both from diagramatics [26] and from symmetry consideration [125] 

L oo 

//f^) = 2/J]J]2Mj)P^,+i, (2.73) 

1=1 j=0 

where "Pf j+i is the projector onto the spin —1 — j representation appearing in the tensor 
product of two spin —1/2 representations and the harmonic sum h{j) is defined by 

^ 1 

^(^■) = E ^ = ^(^' + 1) - V'(l) • (2-74) 

k=l 

where ijj{x) = d/dxlnr{x). This is the s = —1/2 case of the integrable system con- 
structed in [137] and explained in [133]. The Hamiltonian (2.73) is diagonalized by the 
following Bethe ansatz 

Uk + i/2J ^^Uk-ui-i 

Finally, the ]9sm(2,2|4) Hamiltonian at one loop can be diagonalized by the following 
Bethe ansatz 

— (2.76 

M2 .Ml . , M3 

^ ^ TT U2,k - U2,l - I -l-T U2,k - Ui^l + V2 -|-|- 

f4 - '^2,1 + i fj|- U2,k - Ui,l - i/2 f \ 

^ _ TT «3^^^2^W4^^J^_V2 -TT tt3,fc - U2,l + l/2 

f=T ^3,fc - uu + i/2 Ms^fc - U2,i -i/2 

M4,fc + V 2 \ _ -rr tt4,fc - IMd + i TT ^4,fc - ^3J - V2 -|-r 

U4,k -i/2) U4,k - U4,i - i M4,fe - + i/2 ^1 

^ _ J-j- ^5,fc - U4,l - 1/2 j^U5^k- U6,l + l/2 



jJl U5,k - U4,i + i/2 U5,k - U6,i -i/2 

Me -Mr , . M5 

^ ^ TT ^6,fc - Uq^I - ^ TT UQ,k - U7J + l/2 

J^k "^''^ ~ ^ ^ f=l '^^'^ ~ ~ tl 

Me . ,„ 

^ _ TT UT^k - Uq^I + V 2 

r=T - M6,z - i/2 ' 
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U2,k 


- us^i + i/2 


U2,k 


- U3,i - i/2 


UA,k 


- M5,i - i/2 


UA,k 


- U5^i + i/2 


U&,k 


-U5^i + i/2 


U6,k 


- U5,i -i/2 
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Figure 3: Duality transformation on a piece of Dynkin diagram . 
with energy 

M4 

In order to avoid Bethe roots at infinity, the magnon numbers obey Mi < M2 < M3 < 
M4 > M5 > Mq > M7. The equations (2.76) arc written for the Dynkin diagram from 
Fig. lb. For the Lie superalgebras, since the Dynkin diagram is not unique, to different 
Dynkin diagram correspond different sets of nested Bethe ansatz equations. Of course, 
the different sets of Bethe ansatz equations should be equivalent to each other. This can 
be checked on the sZ(2|l) case on the Bethe equations for the supersymmetric t-J model 
[138, 139] by performing the so-called duality transformations on the Bethe roots. The 
general case for an arbitrary root system of su{N\M) was treated by Tsuboi [140]. A 
duality transformation act on a sequence of two nodes as shown in the figure 3. Let us 
exemplify it first for the case Mi = M2 = M5 = Mq = M7 = 0, with equations 

1 = n ''''' ~ ~ ^ /o (2-78) 

, ■ \ L Mi , ■ Ms . ,„ 

U4,k + V2 \ _ -p-r U4,k - + 1 -TT tt4,fc - ^3,; - 1/2 

U4,k U4,k - U4,l - i U4,k - U3^i + i/2 

which correspond to a s/(2|l) subsector and is the case considered in [139]. Let us 
consider the polynomial P{v) of degree M4 — 1 

Mi Ma Ms M4-M3-I 

P{v) ^ Y[{v - U4,i - i/2) - \{{v - U4,i + i/2) = Hiv - u,,,) n - ^3,.) (2-79) 
;=i 1=1 k=i j=i 

As the second equality sign suggests, the magnon rapidities u^^k are among the roots of 
P{v), due to the first set of equations in (2.78). Since M3 < M4 — 1, P{v) may have 
other roots, denoted by v^j. In the same time we have 

Ms M4-M3-I 



P{u4,k + i/2) = -i Yl{u4,k - U4,l + i) = Yl{u4,k - + V2) (^4,/= - ^^3j + 'i/2) 

l^k 1=1 j=l 

Ms Mi-Ms-l 

P{u4,k-i/2) = -iYl{u4,k-U4,i-i) ^Yl{u4^k-U3,i-i/2) JJ {u4,k - Vsj - i/2) 



l^k 1=1 j=l 
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We conclude that the equations (2.78) are equivalent to the set of equations 




V3,j - U4,i - ^2 
V3,j - U4,l + «/2 



U4,k - V3,l + V2 

U4,k - vs^i -i/2 



(2.80) 



In particular, when M4 — 1 = M3, the equations (2.80) become the gl{l\l), purely 
fermionic equations 



which are equivalent with the Bethe ansatz equations for the XY model. Another 
interesting case, which can be treated analogously, is the one with M4 — 1 — M3 = M5. 
Doing a duality transformation on both sides of the central node 4 we obtain 



which are nothing else than the Bethe equations (2.75) for the sl{2) sector. This is a 
consistency check for the one loop Hamiltonian (2.73). 



3.1 The notion of perturbative integrability and the asymptotic 
Bethe ansatz. 

Once that the integrability was established at the level of one loop in the gauge theory 
and at the classical level in string theory, the question arose whether the two integrable 
structures are not manifestations of the same integrable structure, which would exist at 
any order in perturbation theory in the gauge theory, and at the quantum level in string 
theory. The first piece of evidence that this might be the case was put forward by Beisert, 
Kristjansen and Staudacher [27], who computed the two loop order dilatation operator 
in the so(6) sector. Upon diagonalization, the spectrum of anomalous dimensions, in 
the planar limit, showed degeneracies which could not be explained by any symmetry of 
the theory and which are characteristic of integrable hamiltonians. The existence of this 
degeneracy led the authors of [27] to conjecture that the dilatation operator is integrable 
at any loop order. 

In the su{2) sector, the two loop dilatation operator of [27] can be written again in 
terms of the permutation operators as 




(2.81) 




(2.82) 



3 



Perturbative integrability in A/^ = 4 SYM 



L 



L 




(3.1) 



1=1 



i=l 
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and it connects the next-to-nearest neighbors. Prom planarity of the Feynman diagrams, 
we infer that at nth loop order the spin can interact at distance at most n, therefore 
the dilatation operator would be a long range interacting spin chain with the range 
of the interaction controlled by the order in perturbation theory. It is clear that D4 
is not an integrable Hamiltonian by itself, and it does not commute with D2 so the 
two pieces cannot be diagonalized simultaneously. However, one can interpret D4 as a 
piece of the deformation of the integrable Hamiltonian D^. The notion of perturbative 
integrability can be understood as follows. The integrability of D2 implies that there 
are other conserved charges, generated by the expansion of the transfer matrix (2.56). 
Let us call /*^"^ with n = 3, . . . , L the extra conserved charges, in addition to the shift 
operator and the Hamiltonian D2 = 

[j(n)^jM| ^ Q _ (3 2) 

When the coupling constant is turned on, the hamiltonian and the other conserved 
charges get deformed into {D{g) — Do)/g'^ and I'^'^^g) such that 

[I^''\g)J^"'\g)] = 0. (3.3) 

When expanded'^ in powers of g"^, the conserved quantities are not conserved order by 
order in g. A practical way to determine the conserved quantities recursively order by 
order in perturbation theory is to impose 

k k 

E 4fg^^-^ Y: liT'g''-'] = + Oig^'^) . (3.4) 
1=1 1=1 

The authors of [27] used the definition (3.4) and the planar structure of the Feynman 
graph to predict what would be the three loop dilatation operator, in the case this would 
be integrable. Their proposal is 

L L L 

De = -4^(P,,,+3P^+M+2-n^+2P^+l,^+3) + 16 5](Pv+2-l)-56^(P,,,+l-l) . (3.5) 
i=l 1=1 i=l 

It was later confirmed by Beisert [39] that this expression is completely fixed by super- 
symmetry. The three loop dilatation operator was therefore proven to be integrable. 
The most advanced to date perturbative computation which confirms that the dilatation 
operator is integrable up to four loop is due to Beisert, McLoughlin and Roiban [141]. 

Several questions can be asked once we have determined that the dilatation operator 
is integrable beyond the one loop order. One of them is how to diagonalize the higher 
order hamiltonians. The other is how to determine the interaction at orders n > L, when 
the Feynman graphs can be wrapped around the spin chain. 

A partial answer to the first question came [34] from the already known spin chains 
with long range interaction. The spin chain which can reproduce the structure in (3.1) 

^In the su{2) sector, the dilatation operator expands in even powers of g only, because it conserves 
the parity of the number of scalar fields. This should be also the case for the other conserved quantities. 
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and (3.5) is the Inozemtsev spin chain [33]. Unfortunately, the Inozemtsev spin chain 
involves only two-spin interactions and is not general enough to reproduce the many-spin 
interaction which proliferate in the dilatation operator at higher loops. However, using 
the higher order conserved Hamiltonians of the Inozemtsev model, it was possible to 
match the integrable structure of the dilatation operator up to three loop. In particular, 
it was possible to see that the dilatation operator can be still diagonalized by Bethe 
ansatz, provided that the chain is very long. The term asymptotic Bethe ansatz was 
invented to describe this situation. Unlike other integrable spin chains, in this spin chain 
the magnon energies are not additive; the correction to the two-magnon energy vanishes 
exponentially with the length of the chain. 

The comparison with the Inozemtsev model also allowed to check that the Yangian 
symmetry survives at higher loops, at least asymptotically [34]. The Yangian is a sym- 
metry for the Heisenberg spin chain [142] in the limit of infinite length, and an exact 
symmetry for the Haldane-Shastry model. Dolan, Nappi and Witten [143] emphasized 
the link between integrability and the Yangian in the context of AdS/CFT. Agarwal 
and Rajeev [144] checked that the two-loop Yangian coming from the Inozemtsev model 
commutes with the two loop dilatation operator. Beisert and Erkal [145] studied the 
higher loop deformations for the long-range gl{n) spin chains proposed by Beisert and 
Klose [146]. 

Later, we are going to see that the dilatation operator, up to three loops, can also 
be reproduced by a projection of the Hubbard model onto the spin sector. However, to 
date there is no closed expression for the spin chain which corresponds to the spin part 
of the Hubbard model. 

Different long-range integrable deformations of the nearest-neighbor Heisenberg spin 
chaine were considered in the literature, see[146] [147] [148]. 



3.2 The Inozemtsev model 

The Inozemtsev model [33] describes a periodic chain of L spins in the su{2) spin 1/2 
representation which continuously interpolates between the short range Heisenberg model 
and the long range Haldane-Shastry [149, 150] model. The Inozemtsev Hamiltonian is 
defined as 

L L-1 

Hi = J2J2^L,n/.{n){l - P,,,+„) , (3.6) 

j=l n=l 

where the interaction strength V^T^/Kin) is given by the Weierstrass elliptic function with 
periods L and Itt/k, 

jy C2;) — J- -)- f ^ ^ ^ C3 7) 

■^-.tt/^ ^2 — mL — m7r/fi;)2 {mL + imr / k)'^ J ' 

where in the primed sum the terms with m = n = is omitted. When one of the periods 
become infinite, the Weierstrass function degenerates either to a trigonometric or to a 
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hyperbolic function 



limP,,„.W = (£)^(_^-i) . (3.8) 
= + . (3.9) 

The first case corresponds to the Haldane-Shastry model, which played a special role in 
the landscape of integrable models with long range interaction. This is one of the few 
cases where the Yangian is an exact symmetry of the Hamiltonian, even for systems of 
finite length, and it is also an example of system with purely statistical interaction. The 
second case is the one we will be concerned with in the following. 

The Heisenberg model can be obtained in the limit when the imaginary period van- 
ishes, K. ^ oo. We define 

t = e-'^" (3.10) 

such that the interaction strength becomes, after removing an additive and multiplicative 
constant 

Vr / (v) 1 fn+Lm 
KM = - I2 = 2: (,n _ tLmy ■ (3-11) 

The interaction strength (3.11) has the property, characteristic of the Af — A SYM 
dilatation operator, that ht^Lin) ~ In particular, we have 

\im ht,L{n) /t = 5n,i + Sn,L-i (3.12) 

The full solution of the Inozemtsev model is highly involved and the wave functions 
can be found via a correspondence [33] with the elliptic Calogero-Moser system. The 
diagonalization cannot be done via Bethe ansatz and the eigenenergies are not the sum 
of the energies of individual magnons. The situation simplifies when the length of the 
chains is very large and 

t" 1 

hthi'n) = -, TTT = 9 • (3.13) 

In this situation, the solution can be obtained by Bethe ansatz, with 

with the rapidity 0(p) and magnon dispersion relation £(p) 



= n t:: : p-") 



V'W = ^of^V^cfl^) . (3.15) 



27rzK J 22K \2k 

(—\- —V ( ^ 



<P) - 1^CAZ)-T-,'pA'1]-^'\p) (3-16) 
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where Vi{z) — dCi{z)/dz denotes the Weierstrass eUiptic functions of periods 1 and in/K. 
The non-additivity of the spectrum can be compared with that of the Y system (8.60), 
where at finite length there are corrections to the energies of the individual magnons. 

The hyperbolic Inozemtsev model posseses a Yangian symmetry which is a deforma- 
tion of the Yangian symmetry of the (infinite length) Heisenberg chain. The generators of 
the Yangian symmetry can be obtained by analytical continuation of the ones [151, 152] 
of the Haldane-Shastry spin chain, by replacing Zk — e'^'^'^^l^ by z^ — t"^ 

So" - E^"' 

j<k 

with a'^ being the Pauli matrices associated to site j. The existence of a monodromy 
matrix, associated to the Yangian (3.17) allows to compute [151, 153], although relatively 
painfully, the higher order conserved charges of the hyperbolic Hamiltonian 



where the prime indicates that the sum is over distinct summation indices, Zij = Zi — Zj, 
and Pijk, etc. represent cyclic permutations of spins. 

3.3 The three loop dilatation operator and Inozemtsev model 

Since the fourth order conserved charge Hj'^^ of the Inozemtsev model contains a four-spin 
interaction, we can try to use it in order to match the first few orders of the dilatation 
operator 

D{g) = L + h{9)HP{t) + h{9)Hf\t) + ... (3.19) 
Explicit computation shows that if we relate the coupling constant to 

i = /-3(?^ + lV + ... (3.20) 

one can reproduce the terms (3.1) and (3.5). Using the rapidity (3.15) representation in 
terms of the infinite sum 

(^(p) = -cot- + -E ^cot(--iKn)+cot(- + iKn)J (3.21) 

n>0 

2 CO 2 ^ (1 - t^f + 4t« sin2(p/2) ' 



H 



I 



af = J]^^(P,,,-i), 



I 

(4) _ ZiZjZkZi 

I — 2^ [^ijkl 

ijkl ^V^ik^kl^^i 
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one can use the Bethe ansatz (3.14) to compute perturbatively the anomalous dimensions 
up to three loops in the su{2) sector. 

We are going to see later that the same problem can be solved starting from the 
Bethe ansatz of the Hubbard model. Neither the Inozemtsev nor the Hubbard model are 
able to reproduce the fourth order dilatation operator. 

A somehow different approach was taken by Staudacher [38], who proposed to di- 
agonalize the operators (3.1) and (3.5) on an infinite chain by generalizing the proce- 
dure employed in the coordinate Bethe ansatz for the Heisenberg model (2.48) to a case 
where the interaction involves several sites. This procedure, which he named perturbative 
asymtptotic Bethe ansatz allowed to reproduce the scattering phase up to the third loop 
order. This point of view encouraged to concentrate more on the concept of scattering 
matrix and factorized scattering than on the notions of R matrix and monodromy matrix 
which are central to the algebraic Bethe ansatz approach, and allowed to circumvent the 
difficulties related to change of length which occurs once larger sectors, like su{2\3) are 
taken into account. 

In the scattering image, one considers the reference state with a very large number 
of up spins - or operators of type Z in the = 4 SYM case 

\n) = \zzzzz ...zzzzz) (3.22) 

An elementary excitations above this "vacuum" is a magnon, which corresponds to re- 
placing of one of the Z by another field, for example X. One therefore reduces the 
problem of solving completely the problem of L spins to that, simpler, of scattering of 
magnons above the (infinite) sea of Z fields. Generically, the magnons carry quantum 
numbers and their scattering is described by a matrix, while in the particular su{2) 
case the scattering matrix is just a phase factor. For the su{2\3) chain the scattering 
matrix has su{2\2) symmetry, while for the full psw(2,2|4) the scattering matrix has 
sti(2|2) (g) sti(2|2) symmetry. 



4 Strings on AdS^ x 

String theory on AdS^ x can be defined via a sigma model on a coset space. This 
construction is due to Metsaev and Tseythn [20] , following the corresponding construc- 
tion of superstring action in the ten dimensional fiat space by Henneaux and Mezincescu 
[154] . The coset space is given by 

^^t^(2,2|4) 



SO(A, 1) SO(5) ' 

The bosonic part of PSU (2, 2|4) is given by SU (2, 2) ^SU (4), which is locally isomorphic 
to S0{5, 2) (8) S0{6), such that the bosonic part of the coset space is given by AdS^ — 
SO{5,2)/SO{4:,l) times = S0{Q)/S0{5). The coset model provides a natural way 
to couple strings to the Ramond-Ramond fields. 

In order to be able to write down the action, it is necessary to introduce a few basic 
facts. The presentation here follows the general lines of the pedagogical reference [123]. 
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The superalgebra su{2, 2|4) is a non-compact real form of the algebra s/(4|4), which 
can be thought as the algebra of 4|4 x 4|4 supermatrices over the complex numbers, with 
zero supertrace 

(4.2) 

where the 4x4 elements a and b are even and 6 and rj are matrices with odd (grassman- 
nian) elements and STrM = Tra — Trfe = 0. The non-compact real form is specified by 
the relation 

M^H + HM = , (4.3) 

where 

/ I2 \ 

(4.4) 





(h 













-12 







V 





14 



It should be noticed that the identity belongs to the su{2, 2|4) superalgebra, since STrl = 
4 — 4 = and that it is a central element, since it commutes with all the other elements 
of the algebra. Removing the identity reduces the sit(2, 2|4) superalgebra to psu{2, 2|4). 
The latter algebra has no representation in terms of supermatrices. 

The quotient algebra can be characterized with the help of the automorphism gener- 
ated by the supertransposition via 

n{M) = -KM'^R-^ . (4.5) 

The supertranspose M^^ of M is defined by 

M" =(;:;?'). (4.6) 

such that 

(M^r=(_"^ (4.7) 

K = diag((7, a,a,a) , a ^ T V2 = J 0^ ) ' ^"^'^^ 

We deduce that — 1, and on the even subalgebra this automorphism is of order 2. 
One can decompose the superalgebra sZ(4|4) into a direct sum of graded spaces 

g = gW + gW + gi^) + gi^) (4.9) 

where g^'^^ is the subspace of ^ = s/(4|4) corresponding to il{M) = i^M. The fixed 
point g^^^ of the automorphism Vt coincides with the even subalgebra so(4, 1) © so(5). 



and 
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4.1 The Lagrangian 

We consider g an element from the supergroup SU (2, 2|4) and we introduce the associated 
one-form with values in the su{2, 2|4) algebra 

A = = + A^^^ + ^(2) + ^(3) _ (4jQ^ 

This definition insures that is a fiat connection 

d^Ap - dpA^ - Ap\ = Q, (4.11) 

where the indices a, /3 take one of the values a or r. The Lagrangian density is postulated 
to be 

C^-g W^h'^^^TYiA'i^ Af) + «:e°'^STr(A«A5f^)] (4.12) 
with the coupling constant being related to the 't Hooft coupling constant by 

(4.13) 

/i"^ the world-sheet metric and e^'^ = 1. In the conformal gauge h'^^ = diag(— 1, 1). 
Under right multiplication of q with an element f) from 150(4, 1) (g) SO (5), — > 0f), 
all the elements A^'''^ with i = 1,2,3 transform by adjoint action A^"^^ — )■ i)^^A^'^^i), 
while ^4^'^) — >■ \:)~^A^^^ f) — l)~^dl). This means that the string lagrangian is invariant 
under local transformations in 5'0(4, 1) <S) SO{5). Moreover, the lagrangian is invari- 
ant by a shift involving the identity matrix, so it depends on the coset elements from 
PSU{2, 2|4)/.SO(4, 1) ® 50(5). 

The second term in (4.12) can be interpreted as a topological term and it contains the 
fermionic degree of freedom. In addition to the local symmetry generated by the right 
multiplications, the action possesses another symmetry, generated by left multiplication 
with elements associated with the odd part of the psu{2, 2|4) algebra. This symmetry is 
the so-called kappa symmetry and it can be used to eliminate 16 non-physical fermionic 
degrees of freedom and leave only 16 physical fermionic degrees of freedom. The kappa 
symmetry requires the coefficient in front of the Wess-Zumino term to be k = ±1. For 
more details on the kappa symmetry we refer the reader to [123]. 



4.2 Integrability of the classical sigma model 

Classical integrability can be formulated as the compatibility condition for a system of 
partial differential equations in two dimensions (cr, r) 

9,* = L,(a,T,^)* (4.14) 

Here, \1/ is a vector of dimension n and L^r^r are nxn matrices. Both \E' and L depend on 
an extra complex parameter z, the spectral parameter. The system (4.14) is compatible 
provided that the matrices L^-^t- obey the zero curvature condition 

d^Lr - drL„ - [L„, Lr]=0. (4.15) 
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The existence of the Lax connection (4.15) L^r^r allows to define the monodromy matrix 

/•27r 

T{z) = Pexp / daL^{z) (4.16) 

JQ 

where P is path ordering operator and the system is considered to be on a cylinder of 
circumference 2ti. Due to the zero curvature condition, the trace of the monodromy 
matrix is independent of r. 

Construction of the Lax connection 

The classical integrability of the coset model was established by Bena, Polchiski and 
Roiban [25] , who proved the existence of the Lax connection. This is written as an linear 
combination of the different Z4 components of the flat connection 

La = coAW + ciAi^) + c^Ta^e'^Mf + + c4A(f) (4.17) 

Imposing the condition that the connection La is flat is equivalent to the equations of 
motion if the condition for the kappa symmetry is obeyed 

= 1 (4.18) 

and the coefficients q are determined by a single parameter z 

=. = i('~' + ^). =' = -^(-'-^). <=» = ^ = - 

The action of the automorphism Q translates into 

n{L^{z)) - L^iiz) (4.20) 

A different presentation of the spectral parameter, which will be used in the following is 
given by 

z' = ^ (4.21) 

X — 1 

such that 

n{L^{x)) = L,(l/x) (4.22) 

The Lax connection becomes singular at the points x — ±1, or 2; = 0, 00. The mon- 
odromy matrix T{z), defined by the equation (4.16) depends only on the point where 
the contour is cut open. Changing this point will only change T{z) ^{z)T{z)^~^{z), 
such that the eigenvalues of the monodromy matrix encode the physical information. 
Diagonalizing the monodromy matrix we obtain [42] 

u{z)T(z)u-\z) = diag(e'^i(^\ e'^'^'\ e'^'^'\ ^iM^)^^ipi{^)^ ^p2{z) ^ ^p^{z) ^ ^p,{z)^^ ^4 23) 

where the two groups of eigenvalues correspond to the two gradings. Unimodularity 
imposes SDetT(2;) = 1, which translates into the condition 

Pi + P2 + P3 + P4 - Pi - P2 - Ps - P4 e 27rZ (4.24) 
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The eigenvalues e^^^^^\ e'^*^^^ are analytic with respect to z except at the points 2; = 0, oo 
and the points Za where eigenvalues become degenerate. The typical case is that when 
two eigenvalues cross and it can be analyzed on a 2 x 2 matrix 

with eigenvalues 

71,2 = ^(a + d± V(a - + 46c) (4.26) 

The singular point Za corresponds to the point where {a — + Abe — 0, and expanding 
in z around Za one gets 

72 = 71 (1 - OiaV^^a + 0{Z - Za)) . (4.27) 

We conclude that the typical singularity when two eigenvalues of similar grading collide 
is of the type square root. When two eigenvalues of different grading come together, we 
have to analyze the case when F is a 2 x 2 supermatrix, with eigenvalues 

be be 

71 = + a , 72 = + d (4.28) 

a — a a — a 

where b and c are grassmannian variables with (6c) ^ = 0. When a approaches d, the two 
eigenvalues become close to each other. So the typical singularity when two eigenvalues 
of different grading approach to each other is a pole. At z = 0, oo, the functions e'^^^^^\ 
Qtpk{z) j^g^yg essential singularities of the type e"°/^ and e"°°^ . In order to remove this 
singularities, one can consider the logarithmic derivative of the monodromy matrix 

d 

Yiz) = -iz— InTiz) 
az 

whose eigenvalues Uk = zd^ Pk he on an algebraic curve, since their only singularities can 
be poles or branch cuts. 

The automorphism (4.5) acts on the matrix 

y{z) = diag{yi{z),y2iz),ysiz),y4iz)\yi{z),y2iz),ysiz),y4iz)) (4.29) 

by 

^{y{z)) =y{iz) = -dmg{y2{z),yi{z),y4{z),m{z)\m{z),Mz),Mz),y3{z)) (4.30) 

since the conjugation with the matrix a — i~^a2 exchanges the diagonal elements. We 
deduce that 

yk{-z) = ykiz) , yki-z) = yk{z) (4.31) 
so that the algebraic curve can be parametrized by 

^2 + 1 

- = ^ • (4-32) 
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Under this parametrization, the sheets of the algebraic curve are permuted by the trans- 
formation x — )■ l/x as follows 

mix) = mil/x) , yi{x) = mil/x) (4.33) 
m{.x) = y4{l/x) , ^3(2;) = ^4(1/3;) 

It is clear that at the fixed points of the transformation x 1/x, that is x = ±1, the 
algebraic curve has singularities. A particular case where these singularities collide with 
the branch cuts of the type (4.27) will be studied in section 7.3. 

5 Comparison between the gauge theory and the 
string results 

An important point in verifying the hypothesis of integrability at any value of the cou- 
pling constant was the comparison between the results obtained in the weak coupling 
regime, corresponding to perturbative gauge theory, and the strong coupling regime, cor- 
responding to perturbative string theory. A priori, the two regimes do not superpose with 
each other, so it looks like the results cannot be compared. However, Berenstain, Malda- 
cena and Nastase [21] proposed to study a regime where a string moves very fast along 
the equator of 5*5, with angular momentum J. In terms of the spin chain, J — L — M. 
In this regime the effective coupling constant is the so-called BMN coupling constant 

A' = ;^ or g'^^ (5.1) 

If one considers A and J large but A' small, then one may compare the results obtained 
from the regime A small, J large. Of course, nothing insures that there will not be 
order of limits problems, such that the results in the two regimes may differ. Luckily, it 
occurred that the results coincide up to order A'^, and these results were enough to allow 
conjecturing the general structure of the Bethe ansatz equations. Starting with the third 
order in A', the string and gauge theory results are different and at four loop order the 
dilatation operator does not have a well defined BMN hmit anymore. 

5.1 Strings in the plane- wave background. The BMN result 

The plane-wave geometry is the geometry seen by a particle moving at a speed close to 
the speed of light [22] . In the case of the AdS^ x geometry with metric 

dsAdS^xS^ = R^[-dt^ cosh^ p + dp^ + sinh^ p d^l + dijj^ cos^ 9 + dO^ + sin^ Od^'^] (5.2) 

one can consider a massless particle moving on the light-like trajectory parametrized by 

p = 0, ^ = 0, t = t{<f)), il^ = ilj{<f)) (5.3) 

Strictly speaking, the string background is AdSr, x with the otherwise compact coor- 
dinate t decompactified. By introducing the light-cone coordinates = [t ±ip)/2 and 
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rescaling 



= , x~ = fiR^ x~ , r = R p , y = R 9 (5.4) 

A* 

with n a mass parameter necessary to keep the dimensions right. Upon sending the 
parameter R to infinity, the AdS^ x 5*5 metric becomes the pp wave metric [22] 

dsl^ = -Mx+dx- - i^'^if + r^){dx+f + df + df^ (5.5) 

with r and y two four- vectors embedded in AdS^ and 5*5 respectively. The mass parameter 
fjL is related to value of the self-dual five-form by 



+1234 — -C^+5678 



4/1 (5.6) 



which breaks the symmetry of the transverse modes from SO{8) to 5*0(4) x 5*0(4). 
When the five-form field strength vanishes, // — > 0, the plane-wave metric becomes 
the ten dimensional Minkowski metric. The energy E = idt and angular momentum 
J = —id^ get translated, in the light cone variables x^ into 

2p- = id^+^iiiidt + d^)^iiiE-J) (5.7) 
2,- = ,8.-^^0.-8,)^^ (5.8) 

The light cone momentum p"*" is non- vanishing provided J grows in a correlated way with 
R^, which means J ^ a 'VX. The light cone Hamiltonian is given by Hic = 2p . Under 
the AdS/CFT correspondence, the energy E should correspond to the conformal dimen- 
sion A of some local operator in the A/" = 4 SYM theory, while the angular momentum 
J, which is a 50(6) charge, should correspond to one of the R-charges. In particular, 
the case E — J corresponds to a point-like string moving on the equator of the sphere 
S^. This is a BPS state with l-Lic — and its gauge theory counterpart is a state of the 
type Tr Z'^ . 

The gauge freedom for the reparametrization of the string can be fixed by choosing 
the light cone gauge = r, with r the world-sheet time. 

The action of superstrings in the plane- wave background in light- cone gauge was 
written down by Metsaev [155], who showed that it corresponds to eight massive bosons 
and eight massive fermions, whose masses are equal due to supersymmetry. As it can be 
deduced from (5.5), this action reads 



111 

-{dT-zf {drrzf' z"^ + fermions 

2 2 2 



(5.9) 



with the eight- dimensional vector z = {f,y) and £ = 2Txa'p^ [156]. The action (5.9) can 
be quantized similarly to the free bosonic string action, its spectrum being spanned by 
massive bosonic (fermionic) oscillators ol\^ [9]^, 9^) with energy 
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Including the zero modes which arise in the curved background, the energy of the string 
in the hght cone gauge is given by 

El, - N^Eo + ^(iV„ + Nn)Er, . (5.11) 

n>l 

where the occupation numbers A^„, obey the level- matching condition (or Virasoro 
constraint) 

5^(iV, - iV,) = . (5.12) 

n>l 

Taking into account the fact that £■ ~ J, and i?^ = Aa'^, we obtain from (5.8) 

4 ^ A' (5.13) 



such that the dispersion relation of the impurities (5.9) take the so-called BMN form 



j \„2 

Ejfi = Vl + A'n2 = y 1 + _ . (5.14) 

Let us consider the two-impurity case with Nn = Nn = 1 for some mode number n and 
compare the leading result in A' for the energy of the string and the corresponding result 
in gauge theory. Prom (5.7) we conclude that 

A - J = Eijfx = 2Vl + A'n2 = 2 + AV - ^ + ^ - 0(A'2) (5.15) 

4 8 

In the same time, from the Bethe ansatz equations (2.68) one has, for ui — —U2 — u, or 
equivalently pi = —p2 = p, 

^V" = l, or e-<-') = l, (5.16) 
u — 1/2 J 

with the solution 

The one-loop prediction in gauge theory is then 

A . 9 Tin ^ Xn^ A, . 

^ = ^ + ^ l3T = i + (TTTp + «'(^-') 

where the second relation holds for n <^ L. Comparing (5.15) and (5.18) wc conclude 
that ,/ = L — 2 for two impurities and the two expressions agree up to terms of order 
J^^ at one loop. One can go beyond the one-loop order in computing the conformal 
dimension of the two-impurity operator, either by direct diagonalization of (3.1) and 
(3.5) [27], or by using the Inozemtsev asymptotic Bethe ansatz [34]. The result is 

A(A) = J + 2 + A^(^) + A^E^^) + A^E(3) + . . . (5.19) 
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with 



= ^sin^^, (5.20) 
TT^ 2 V4 J + 1 2y ' 

^ = ^'^^tU jTi + (TTIp J- 

which agrees with (5.15) in the hmit n <^ J and at leading order in J. The 1/J 
corrections to quantum strings beyond the plane-wave geometry were computed in [157, 
158]. While the one and two-loop part of these corrections agree with (5.20), the three- 
loop part disagrees. This was the first occurence of a mismatch beween the BMN limit 
in strings and gauge theory, and it was finally understood to be a consequence of the 
non-commuting sequences of limits which are considered. 

The string solutions considered above are a very particular case, since just one of 
the angular momenta is large, and the string excitations are small fluctuations around 
the pointlike string solution. More general string solutions were considered by Frolov 
and Tseytlin [73, 35, 159, 36], where strings can have large angular momentum in some 
planes in or/and in AdS^. Some of these solutions still depend on the BMN coupling 
constant A', but now the approximation of dilute impurities which is implicit in the 
BMN computation is not valid anymore. The first comparison between the spinning 
string solutions and the one-loop Bethe ansatz equations was done in [160] and it was 
pursued in an important number of papers [43, 44, 45, 46, 47, 48, 49]. The 1/J corrections 
were also extensively checked for these solution. 

An alternative way to compare the predictions of the string theory with the results 
from the spin chain was supplied by the formulation of the (finite gap) solutions of the 
string sigma model in terms of algebraic curves [40]. The continuum limit of the Bethe 
ansatz equations gives also rise to an algebraic curve, very similar to that of the sigma 
model. The comparison between the two curves was very fruitful for fixing the general 
structure of the all loop Bethe ansatz. In the next section, this comparison is presented 
for the particular case of the su{2) sector. 

5.2 Algebraic curves and the continuum limit of the Bethe 
ansatz equations: the su{2) principal chiral model 

The full string sigma model (4.12) is relatively complicated to study in full generality. 

At the classical level, one can reduce the dimensionality of the space in which the string 
moves by a restriction which is similar to the restriction to sectors in perturbative gauge 
theory. We are then left with a (string) sigma model on a subspace of the original 
coset space. While this reduction is classically consistent, it is not consistent quantum 
mechanically. The string, although it moves only inside the restricted space, fluctuates 
on all directions of the original space. However, we can forget for the moment about 
the quantum corrections and study strings which move, for example, only on the sphere 
C via the x R sigma model [40]. Let us consider the coordinates $i, $4 with 
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^ • $^ = 1 be the coordinates of the sphere 5"^, which is the also the group manifold for 
SU (2) , and $o the timelike coordinate 



0H l'Vt )^i% i ]^SUi^ (5.21) 



-$3 + i$4 $1 - i$2 J V ~^ ^ 

The string action is given by 



Ss^ = g j\a j dT[{d^^if-{d^%f] (5.22) 



27r 



—gj da j dr 







^Tr(s-'a„S)^ + (a„$o)' 



(5.23) 



The global symmetry of the action is given by SU{2)l^SU{2)r ~ 50(4). The two SU{2) 
symmetries correspond to the left and right multiplication respectively, with currents 
given by la = daQQ~^ ^^^^ 3a = d~^daQ- Under the first group, {Z, —X) and {X, Z) 
transform as doublets, while under the second {Z,X) and {—X, Z) arc doublets. As the 
notation suggests, the two SU{2) groups correspond to the ones in perturbative gauge 
theory which permute the complex scalars, generated by combinations of the oscillators 
c^, and c^, respectively. In the normalization where the charge of the states in the 
doublets arc ±1, an operator of the type Ti Z^~'^ X'^ + . . . will have a charge L under 
the SU{2)l and charge L-2J under SU{2)r. 

This time we consider the action in the gauge $o = and the energy of the string 
is given by 

poo 

A = 2g dadr^o = Airgn . (5.24) 



The equations of movement for the action (5.22) can be written, in the light cone coor- 
dinates cr± — I (r ± cr) , d± = 8-^ ± 

9+i_-9_i+ + [i+,i_] = 0, (5.25) 

where the first equation corresponds to the current conservation and the second line 
is the equation of motion. These equations are to be supplemented with the Virasoro 
condition, which express the vanishing of the stress-energy tensor, 

llV jl = ^Tr £ = -k' . (5.26) 

The equations (5.25) can be reformulated as the zero curvature condition [161] 

9+ J_ - 9_ J+ + [J+, J_] = (5.27) 

for the connection 

J, = ^ . (5.28) 
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Of course, this connection is a particular case of the one presented in section (4.2). The 
monodromy matrix 

QW = Pexpj^"^(^-^) (5.29) 

can be diagonahzed in terms of the quasi- momentum p{x) defined by 

TrQ(x) =2cosp(x) . (5.30) 
Near the singular points x = ±1, the quasi-momentum behaves as 

pW = --^ + .... (5.31) 

while the asymptotics at a; = 0, oo are related to the conserved charges, see [40] for more 
details, 

p[x) = -- — — + ... oo (5.32) 

2gx 

L 

p(x) — 27rm + —X + . . . x ^ . 

This behavior suggests that the transformation x ^ 1/x exchanges the two su{2) com- 
ponents. In order to remove the poles at x = ±1, one introduces the resolvent 

G(x)=p(x) + -^ + -^. (5.33) 

which is analytic on the physical sheet, with possible branch cuts. The resolvent can be 
represented in terms of a density p{x) with support on the branch cuts of G{x) 

G{x) - jdy^y (5-34) 

such that 2'7iip{x) = G{x — iO) — G{x + iO). The asymptotic conditions for p{x) at x = 
and oo translate into normalization conditions on p{x) 

dx p{x) — — , (5.35) 



x^ 2g 
The continuous part of p{x) on the cut is 

p{x + iO) + p{x - iO) = 27mc (5.36) 
where nc is specific to the branch cut C, or, in terms of the resolvent 

G(x + iO) + G(x -z0)=2-f dy^^ = 4^ + 27rnc • (5.37) 

J X — y x^ — 1 
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This equation, together with the normahzation conditions (5.35) can be compared with 
the continuum hmit of the Bethe ansatz equations (2.61). Let us take the logarithm of 
(2.61) and rescale the rapidities Uk — > Luk 

fl + i/2Luk\ ^ ^ f l + i/L{uk-Uj)\ 
PkL = L]n[- TT^Ti— =27rznc+> In ^ ^ 5.38 

For L large, expanding in 1/L one gets 



- = 2nnc + 2 f dv^-^ 
u J u — V 



(5.39) 



with f du p{u) = J/L. In order to compare with (5.37), we make the redefinition 
X — >■ xL/g and take into account (5.24) that At^k = A/gr such that (5.37) becomes 



= ^^^ + 2.n.. (5.40) 



This transformation introduces explicitly the BMN couphng constant 

A' = leTT^^VL^ = IGtt^^'^ . 

Since A = L + O(A'), we see that at leading order in A' (5.40) and (5.39) coincide. We 
can go with the comparison at higher loops, by using the Inozemtsev asymptotic Bethe 
ansatz (3.14). If we replace the Inozemtsev rapidity tp Lu, we obtain perturbatively 
the potential on the l.h.s. of (5.39) 



pL = - + 2g'^^ + 6g'^^ + ... (5.41) 



dx[l-^—] p{x) = -= I du p{u) (5.42) 



In this expansion, we dropped the higher powers of 1/L. The terms under the dots do not 
have the so-called BMN scaling, in powers of gi' = g/L plus 1/L corrections, anymore. 
On the other hand, from the first and the last normalization conditions (5.35) we deduce 
that 

the last equality being the normalization condition for the density in the spin chain. The 
relation (5.42) suggests that the Bethe ansatz rapidity u and sigma model rapiditiy x 
are related by 

u = X ^ , pix) = piu) (5.43) 

X 

at least perturbatively. Inserting this definition in the perturbative expansion (5.41) of 
the potential we obtain 

pL=- + g''-^+g'^'\ + ...^^^^. (5.44) 
X x-^ x° x^ — g'^ 
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The proposals (5.43) and (5.41) were made by Beisert, Dippel and Staudacher (BDS) 
[50] in order to extend the Inozemtsev Bethe ansatz beyond the three loop order. The 
corresponding equations, written in the variable are 

J x-y x^ - g'"^ x^ J (x - y)(l - g'^/xy) 
while the corresponding equation obtained from the sigma model reads 

2 / dy^^^ = ^ + — / dy ^'^^^ + 27rnc (5 46) 

J x — y x'^ — g''^ x"^ J {x — y){l — g''^/x'^) 

The two equations are almost identical, except for a term in the denominator in the 
integral in the r.h.s. Since the two equations are different, we have assigned different 
indices to the corresponding densities. When expanded in powers of g', they agree up to 
order g'^ and start to disagree from order g"^. The mismatch will occur in the conformal 
dimensions at order g'^. At order g^, both Inozemtsev and BDS Bethe ansatz are known 
to be correct, therefore the mismatch comes from a different order of limits taken to 
arrive at the two formulas. In the case of the string sigma model (5.46), one takes 
g — )■ oo, L ^ oo with g' = g/L fixed, then g' 1 while in the case of the gauge theory 
spin chain (5.46) one takes g finite, L ^ oo and then g' = g/L ^ 1. Knowing the 
all-loop Bethe ansatz allows to remedy this mismatch. 

The comparison between the BDS chain and the su{2) principal chiral model pre- 
sented above involves only one of the su(2) symmetry sectors. A spin model with two 
su{2) types of spin was proposed by Faddeev and Reshetikhin [162] as a discretization 
of the principal chiral model, see also [133] . 

The algebraic curve for the whole system was obtained in [42] . Intermediate steps in 
obtaining this result consisted in analyzing the thermodynamical limit of the so{6) [163] 
and psw(2, 2|4) [164] one-loop Bethe ansatze. 



6 All loop integrability in A/" = 4 SYM 

The first step in extending the Bethe ansatz to all loop order was made by Beisert, 
Dippel and Staudacher [50], who proposed how to extend the three-loop su{2) Bethe 
ansatz obtained from the Inozemsev model such that it matches at least partially with 
the prediction of the string sigma model and with the BMN dispersion relation (5.14). 
As it was explained in the previous chapter, the BDS ansatz is not able to reproduce the 
results at strong coupling. There is therefore a missing piece which allows to interpolate 
between weak and strong coupling, and this piece was called dressing phase, for reasons 
which will be clear later. It was by no means obvious that the all-loop ansatz proposed 
by BDS is consistent and that it corresponds to any spin chain, be it with multiple 
spin, long range interaction, except from the fact that it matched a fifth-order integrable 
extension of the dilatation operator. 

Shortly after, Beisert and Staudacher [52] conjectured the sm(1|1) symmetric ^'-matrix 
for the sm(2|1) sector at all loops, and, from the consistency relations which the nested 
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Bethe ansatz has to satisfy, tlicy conjectured the full psu{2, 2|4) Bethe ansatz at all loops. 
As in the case of Inozemtsev spin chain, the Bethe ansatz is only asymptotic, that is, is 
valid for L ^ oo. 

A new step forward was made by Beisert [53], who analyzed the s-u(2|2) symmetric 
^'-matrix of the su{3\2) spin chain. This sector corresponds to one of the two branches 
of the Dynkin diagram and is also the smallest sector in which the spin chain can change 
its length. While on a finite chain this property seems to introduce an insurmontable 
difficulty, on a chain of infinite length it allows to modify the magnon symmetry group 
from sm(2|2) to centrally extended s-u(2|2), the extra central charge being related to the 
magnon momentum. This symmetry allows to determine the magnon dispersion relation, 
which is identical to that proposed by BDS [50]. The same symmetry imposed to the 
(S'-matrix, which is supposed to be factorizable, fixes the (S'-matrix up to an overall factor 
called the dressing factor. As explained in chapter (2.3), the choice of the reference 
state breaks the psu{2, 2|4) symmetry to s-u(2|2) (8> s-u(2|2), with the two factors playing 
symmetric roles. Therefore, imposing the symmetry is enough to find the full ^'-matrix, 
modulo the dressing factor. The Bethe ansatz associated with the (S'-matrix determined 
by Beisert [53] coincides with the one conjectured by Beisert and Staudacher [52]. The 
only assumptions arc that the system is integrable and the asymptotic condition L — )■ oo. 

Although the A/" = 4 SYM Bethe ansatz equations where conjectured and then proven 
using symmetry considerations, a concrete realization of the corresponding spin chain 
was lacking. It was surprising to find that the BDS ansatz can be obtained by projecting 
the Hubbard model at half filling, in the way in which the Heisenberg model can be 
obtained from the Hubbard model at half filling in the limit of strong Coulomb interaction 
[165]. The one dimensional Hubbard model is integrable and its solution in terms of 
Bethe ansatz was found by Lieb and Wu [166]. In the language of the Hubbard model, 
the — A SYM magnons are bound states of fermions. In this representation, it is 
straightforward to obtain the BMN/BDS dispersion relation for the magnons. 

The relation between the Hubbard model and the A/" = 4 SYM magnons is deeper, 
and it fair to say that it is not yet completely understood. It was noticed by Staudacher 
[167] and elaborated by Beisert [55] that the 5*- matrix with centrally extended su{2\2) 
symmetry is essentially the same as the i?-matrix of the Hubbard model [168] . 

6.1 The BDS conjecture 

The computation the dilatation operator at higher loops was done perturbatively up to 
five loops [125] in the su{2) sector. The constraint of integrability, the structure inherited 
from the perturbative gauge theory and the condition of BMN scaling fixed uniquely the 
Hamiltonian of the spin chain. However, it became more and more clear^ that the BMN 
scaling should not hold at arbitrary loop order. If one relaxes the BMN constraint, more 
integrable spin chains exist. 

Let us concentrate first to the chain which obeys the BMN scaling at five loops. 
Beisert, Dippel and Staudacher [50] conjectured that this spin chain exists at any loop 

^See the discussion in the previous chapter about the three loop mismatch. 
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order and it is diagonalized by the following Bethe ansatz 



--=TT 44^^44^ (6.1) 

J- A 111 Hi.) 7(1 Hi) — i 



:^ u{Pk) - u{pi) - i 



with 



and the energy given by 



K(p)--cot|Wl + 1652sin^| (6.2) 



^^^s = E (^i + ieTsi^ - (6.3) 

The equations (6.1) should be supplemented with the condition of translational invari- 
ance 

M 

Yl e'P>' = 1 . 

k=l 

As wc have seen in the previous chapter, one can introduce another rapidity variable x, 
related to u by the relation^ 

u^x + — (6.4) 

X 

In this notation, the momentum p can be expressed as 

r+ 

e'f = — , with x^ = x(u ±i/2) , (6.5) 

X 

and the BDS equations can be rewritten as 

+ \L M 

^"1=17 "^h^llll (6.6) 
with the energy 



The higher conserved charges of the BDS spin chain, which generalize the higher con- 
served charges of the Heisenberg model, can be put in the form 



M 

= X] — r 

^-^ r — 1 

k=l 



(6.8) 



The limit r ^ 1 gives the momentum p and r = 2 is, up to a constant, the energy Ebds- 



^In order to avoid too cumbersome notations, we keep the same notation for the rapidities u and x 
in different normahzations. It should be clear from the context which normalization is used. 
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In addition to reproducing the known spectrum of the dilatation operator up to three 
loops, the BDS ansatz has also the merit to reproduce the BMN dispersion relation (5.14). 
For two magnons with pi — —p2 = p ~ 2Trn/L with n <^ L, the energy (6.3) is given by 

Ebmn = 2(yi + A'n2 - l) . (6.9) 

The dispersion relation fixes the two-spin part of the interaction, or more precisely it 
fixes the part of the interaction which is linear in the operators 

Qi,i+n — Pi,i+n 1 

as it was shown by Ryzhov and Tseytlin [169]. The higher power in non-overlapping 
permutations (5i,i+n vanish on the one-magnon states. The corresponding interaction 
strength is simply the Fourier transform of the magnon dispersion relation (6.3) and it 

reads 

Tin 1/2) ^ / 11 



h{n, g) = (^n - -, n + -, 2n + 1, -16/ j (6.10) 

As a function of the coupling constant g, the interaction strength has singularities at the 
points (7 = 0, oo and ±i/4. At these last two points, the interaction strength simplifies, 
to give 

M".±'/4) = ti):_2_ (6.11) 

This is an antiferromagnetic interaction, as opposed to the ferromagnetic interaction we 
retrieve at > 0. This property is also confirmed by the magnon dispersion relation 
which become at this point 

£;(p)|,=±,/, = |cosp/2|-l (6.12) 

which has the minimum at p = tt, which is the antiferromagnetic point. Let us also notice 
that the magnon energies (6.12) are negative, a signal of instability of the ferromagnetic 
state, which will not be the ground state anymore. 

Due to the transformation properties of the hypergeometric function, the strong 
couphng limit of the interaction (6.10) is similar to that at g'^ — —1/16, except from 
the sign of the interaction 

M^,/|.^oc = -^^ (6.13) 
vrn^ — 1/4 

while the magnon energy becomes, at least for l/t?, 

^(P)|.^oo = 4^|sinp/2|-1 (6.14) 

It is interesting that the interaction (6.11) and (6.13) ressembles the inverse-square 
interaction which is familiar from the study of the spin chains like Haldane-Shastry 
[149, 150]. The expressions (6.11) and (6.13) hold for L oo, while for a finite chain 
they should be periodized, such that, for example 

g siv? J 

h{n, 9, L)^^^^ = - . .(2n+i) ■ ^ (5^ (6-15) 
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This interaction differs from the Haldane-Shastry interaction only by the point-sphtting 

. 2 Trn . 7r(2n+ 1) . 7i(2n - 1) 

which is responsible for changing the dispersion relation from quadratic to trigonometric. 

The comments above hold for any spin chain with the dispersion relation (6.3), in- 
tegrable or not and diagonahzed by the BDS equations (6.1) or not. In order to be 
integrable, a spin chain of the type (6.11) or (6.13) has to be supplemented with higher 
spin interaction. It is interesting to know whether such a model could be obtained from 
an integrable finite-difference Hamiltonian like Ruijsenaars- Schneider [170] in the same 
way the Haldane-Shastry spin chain can be obtained [171] from the Calogero- Sutherland 
model. Another common point between the BDS model at strong coupling and the 
Haldane-Shastry spin chain is given by the phase shift 

or N 1, u{pi) - u{p2) + i g-^ . , s , . 

0{pi,P2) = T In -Y—. : — > TT sign (pi - p2 6.16 

I u{pi) - u{p2) - I 

where we considered p G [0, 27r] and considered the determination of the logarithm such 
that the phase shift is antisymmetric, 9{pi,p2) = —0{p2,Pi)- 

6.2 BDS spin chain and the Hubbard model 

A computation which can be easily done for the BDS system in the thermodynamic 
limit L — )■ oo is that of the energy for the antiferromagnetic state, the state with the 
maximum number of magnons M — L/2. This is the equivalent of the problem solved 
by Hulthen [172] in 1938 for the Heisenberg model and amounts to solving the integral 
equation for the density of magnons 

2. +2 fd.- "'"2 , (6.17) 



du {u — u'Y -\- 1 ' 

with 

dpiu) . d , x~^(u) i i 

-«^log — K = ^= , (6.18 



du du x-{u) - Ag^ yju^_ - Ag^ ' 

The solution to the equation (6.17) can be found by Fourier transform and it reads 

P(«)= r (6.19) 

Jo 27r cosh (I) 
wich gives for the energy of the antiferromagnetic state 



EAF{g) = 2gL I — ^-p-^ . (6.20) 



Up to a sign and an overall constant, this is the energy of the antiferromagnetic ground 
state of the Hubbard model at half fiUing computed by Lieb and Wu [166]. This suggests 
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that there is a hnk between the BDS spin chain and Hubbard model at couphng U = 1/g 
and t — —g. 

The Hubbard model [173] is a lattice model for itinerant fermions with spin 1/2. The 
Coulomb interaction is modelled by an onsite repulsion. In one dimension we have 

L L 

^Hubbard = XI + 4+1,^^,^) + ^ ^ J] A-\^iA,\.^i,i " (6-21) 

i=l cr=t,4- '=1 

The operators and Cj^o- are canonical Fermi operators satisfying the anticommutation 
relations 

{ci,a.c,,.'} = {4,sV} = 0, (6.22) 

Due to the Pauli principle, there are four states by site: with no fermions (o), with one 
fermion with spin up (f), with one fermion with spin down (J,) and with two fermions of 
different spins (|). Therefore, the Hilbert space of the Hubbard model is of dimension 
4^, much larger than the Hilbert space of the spin chain, of dimension 2^. A way to 
separate the spin states t, i, from the others is to go at infinite Coulomb interaction and 
half-filling, iVfcimions = L. In this limit, the states with exactly one fermion of either spin 
per site have finite energy, while the states where there is at least one double occupancy ^ 
have infinite energy. Following these states by continuity in the perturbation parameter 
g — 1/U will allow to give an all-loop definition to the BDS spin chain. 

Explicit computation show that the projection of the Hubbard hamiltonian (6.21) 
coincides with the BDS spin chain only up to order g^ on chains with even length. 
This is in contrast with the expectation that the BDS ansatz is correct up order g"^^. 
This slight mismatch can be corrected by introducing twisted boundary conditions for 
the fermion operators or, equivalently, introducing an Aharonov-Bohm flux through the 
system, 

L L 

(e^^^ 4.^^+1,- + ^'"^^ 4+i,.Q,.) - 4tc^t4;ca , (6-23) 

1=1 o-=t4 »=i 
with the twist given by 

= CT^ti (6.24) 

The overall sign of H is opposite to the sign of the usual repulsive Hubbard model; in 
our case we are interested in getting a ferromagnetic ground state, at least in the spin 
sector, 

= in ■ ■ ■ tt) = ci^4t ■ ■ • 4-it4t lo) (6-25) 

We ignore for the moment the fact that the doubly occupied states have, for this choice of 
sign of the Hamiltonian, lower energy than the singly occupied ones. What is important 
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for us is that, at least for ^ 0, the spin states are clearly separated from the states 
with double occupancy. 

The Hamiltonian (6.23) is covariant under particle-hole transformations on the fermi- 
onic operators of either spin [174]. Such a transformation is particularly convenient for 
states close to the BPS ground state (6.25), where it largely diminish the number of 
fermions. We therefore choose to transform 

o ^ t (6.26) 
; ^ I (6.27) 

or, explicitly in terms of creation/annihilaton operators 

Ci,o = 4^ , 4o = Ci,t > (6-28) 

Ci,t = Q4 , = . (6.29) 

This transformation changes the grading of the states transforming a state with an odd 
number of fermions into a state with an even number of fermions, so we may call it a 
duality transformation in the sense of section (2.7). In view of the link, not yet fully 
clarified, of the Hubbard model with a S'u(2|2) symmetric model [55], this transformation 
can be viewed as a supersymmetry transformation. The dual Hamiltonian is 

L L 
H = 9Y.Y1 (^^^^ 4.C^+l><^ + ^'"^^ 4+l,aC^,a) " ^^(1 " 4oQ,o)4^Q,j . (6.30) 

1=1 (7=0,1 1=1 

where 4>^ = (p^, while 0o = tt — (p^. Comparing the two expressions (6.23) and (6.30) we 
conclude that under the duality transformation, the Hamiltonian (6.23) transforms as 

H{g; 0, 0) ^ -H{-g; vr - 0, 0) - M (6.31) 

where M is the initial number of fermions with down spins. 

The one dimensional Hubbard model was solved by Lieb and Wu [166], who showed 
that the coordinate Bethe ansatz was essentialy the same as that for the two component 
fermion gas with delta interaction, problem which in turn was solved by Yang [175] and 
Gaudin [176]. The Lieb-Wu equations for the Hamiltonian in the form (6.23), specialized 
at half-filling, are 

jJ^Uj-2gsm{qn + (f))+i/2' ^ ^ 

'u,-2gsinig^^ + <P)+^/2^^u,~u,+^^ ;t = l,...,M (6.33) 

n=l '"'^ ~ ^^"^{Qn + (p)-l/2 Uk - Uj - 1 

3+^ 



where qk are the momenta of the L fermions and Uj are the rapidities of the M reversed 
spins. The energy is given by 



E^2gJ2cos{qk + (f)) . (6.34) 



k=l 
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The dual Hamiltonian can be diagonalized, by the dual Lieb-Wu equations, which in the 
particular case of half fiUing read 



' -lU,-2,sin(,„-0)+V2' ^-l'---'2^ (6.35) 

'^ u,-2g sin(g. - 0) + i/2 ^ i _ ^g^g^ 

Z\ '^k - 29sm{qn - (p) - 1/2 Uk - Uj - i 

The number of fermions was reduced from L to M + L — {L — M) — 2M and the energy 
is now given by 

2M 

E^-2gY^ cos(gn - 0) - M . (6.37) 

n=l 

The minus factor after the equality sign in (6.36) comes from the particular value of the 
twist. The Lieb-Wu equations with arbitrary value of the twist are written in appendix 
C of [54], using a result by Yue and Deguchi [177]. 

The equations (6.35) and (6.36) allow bound states, which were classified by Taka- 
hashi [178]. Here, we are interested by the particular case when a pair of fermions and 
one magnon make a bound state. Let us first consider the case when just one magnon is 
present and put 

gi-0 = i(7r+p + i/3) , q2-<P=^{7T + p-tP) . (6.38) 

with real p and ^ > 0. Because of the imaginary part of the momentum, at large L the 
l.h.s. of (6.35) will be zero or infinite for n — 1,2. This implies that the r.h.s. has a zero 
or a pole, respectively, 

z z 
u - -^2g sin(gi - 0) , u + -=2g sin(g2 - 0) ■ (6.39) 

This implies that the real and imaginary part of the momentum are related by 

sin - sinh — = — (6.40) 
2 2 4c/ ^ ' 

and we obtain for the magnon rapidity 



«=2Cot|yi + 16y2sin2^. (6.41) 

The magnon energy is given by the sum of the two fermion energies and it also takes the 
BDS form 



E^2g{ sin + sin j - l = Jl + 16^^ ^^^2 | _ ^ (g 43) 
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The bound states are described by the position of their centers u (6.41). When several 
such bound states are present, Lieb-Wu equations can be reduced to equations for their 
centers Uk- The first Lieb-Wu equation becomes 

M . 
fJl Uk -Uj-l 

which is nothing else than the BDS equation, while the second Lieb-Wu equation is 
identically satisfied. The finite size correction are of the order e~^^. At — > 0, we get 
from (6.40) that ^ ~ — Slng^ such that the finite size corrections are of the order g^'^, as 
expected. 

In conclusion, the connection with the Hubbard model provides two key ingredients 
of the all-loop Bethe ansatz, which are the rapidity-momentum relation (6.41) and the 
magnon dispersion relation (6.42). 

It was first noticed by Janik [61] that the variables defined in (6.4) and (6.5) live 
on a torus. A related parametrization, obtained by a Gauss-Landen transformation from 
that of [61], appeared independently in [55] [179] and [78] 

/ V 1 cn s , , 

with the modulus k of the torus 

= (6.45) 

1 + 16^2 V ) 

The torus degenerates at the points k'^ — 0, 1, oo, which correspond to the singular 
values — Q^oo and —1/16 we already encountered in the previous section. Under this 
parametrization, the momentum and the size of the bound states (6.38) correspond to 
the elliptic amplitude function 

p(s) = TT - 2am(K - s,A;), /3(s) = -ztt - 2i am(iK' - s. A;). (6.46) 

and they are related by a translation of quarter of period in both directions on the torus 

p(s + K-iK') = 27r-i/3(s). (6.47) 

6.3 The dressing phase 

As explained in section (5.2), although the BDS proposal reproduces the spectrum of 
the dilatation operator up to three loops in the su{2) sector, it does not reproduce the 
strong-coupling limit obtained from the su{2) principal chiral model, c/. (5.45) and 
(5.46). For the BDS ansatz to reproduce the strong couphng behaviour (5.46), it has to 
be supplemented with an extra piece ^^(xfe, xi) which was called dressing factor in [52] 

f:^V= fi ^^^^^-^(-.,-0 (6.48) 
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At weak coupling, the dressing factor should be equal to one at least to the third loop 
order 

a\xk,xi)^l + 0{g'') (6.49) 

Arutynov, Frolov and Staudacher [51] showed that , in order to reproduce the leading 
strong coupling behaviour, one can assume 

oo 

a'^{xk,xi) = exp2i^g'^''^'^(^qr+2{xk)qr+3{xi) - qr+3{xk)qr+2{xi)^ , (6.50) 

r=0 

where qr{x) are the conserved charges defined in (6.8). This cannot be the right answer 
at weak coupling, since the first corrections at weak coupling appear at order and 
not g^ as they should. The answer to the problem of interpolation of the dressing phase 
between weak and strong will be given in section 6.5. 



6.4 The all-loop asymptotic Bethe ansatz 
6.4.1 The conjecture 

One of the key steps to extend the BDS ansatz to the full psu{2,2\4) algebra was to 
write the equivalent ansatz for the sm(1|1) sector, with spin states Z and ip. Beisert and 
Staudacher [52] verified that the three loop dilatation operator computed by Beisert [39] 
could be reproduced with the ansatz 

f^lV^fl \- a\x,^x,) (6.51) 

Embedding the su{l\l) sector into a larger su{2\l) sector, with spin states Z, X and ■0, al- 
lowed to check that corresponding, su(l|l) symmetric, S matrix satisfies the Yang-Baxter 
equation and therefore corresponds to some consistent integrable model. Imposing the 
other consistency conditions coming from thepstt(2, 2|4) algebra, Beisert and Staudacher 
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arrived at the following ansatz, which generalizes the one- loop Bethe ansatz (2.76) 



1 = 



M2 , M4 

n Ul,k - U2,l + ^/2 -|-|- 1 



U2A 



M2 . Ml 

n W2,fc - U2,l - I -TT U2,k 
l/r. . l/r. 7 -I- ■? J- J. 



V2 f^j 1 - g^/xi,kX4^i 
ui,i + i/2 



(6.52) 



Mi 



U2± 



M3 

n 



i/2 fJ- W2,fc - u^,i - i/2 



-j-j- a^3,fc - a:^; M3,fc - U2,i + i/2 



^ U3,k - U2,l - V2 



1 = 




' M4 

n 



M4,fc - U4,l + i 2 



-n 



1=1 ^4,fe 



•"3^5, fc 



M4 
Me 



U5,k - U6,l + V2 



n 



2^5 



n 



1 = 



1 



n 



Mr 



«6,fc 



n 



M6,« - V2 

+ V2 



«6,z + ^ f-^ iKk 



Mr, 

n 



U(i,k - + V2 

U7,l - i/2 M6,fe - U5,l - V2 



with conserved charges 



n U7,k - U6,l + V 2 T-r ^~ 9 IX7,kXi^i 
, _ W2 11 1 



^/2 fjj; 1 - 9^l^i,kX^^i 



M4 

E 

fc=i 



4,fe; 



4,fc^ 



^r-l 



(6.53) 



These equations should be again supplemented with the zero momentum condition 

M4 



fe=l ^4,fe 



(6.54) 



It is interesting to note that the connectivity of the roots of the Dynkin diagram is 
enhanced compared to the one-loop ansatz (2.76), since the roots number 1 and number 
7 are connected to the root number 4 by a factor which goes to unity when the coupling 
constant vanishes, fif ^ 0, see figure 4. The symmetry of the diagram is also enhanced, 
as the nodes 1 and 3 and 7 and 5 can be interchanged by transforming the x\ jx^ 
and X7 <H- g^ jx^^ at the expense of exchanging L ^ L — Mi + Afs and L ^ L — My + Afs 
respectively''. In fact, as there is no self-interaction for the odd nodes, the roots can be 
transformed one by one from type 1 to type 3 and from type 7 to type 5, Xij g"^ /x^^i 
and X'jj — > g^ jx^^i. Such an operation changes the length of the spin chain L ^ L — 1. 



''Extra factors of — 1 appear for M4 odd. 
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Figure 4: Connectivity of the nodes of the Dynkin diagram for the all-loop Bethe ansatz. 
The dotted hnes represent the nodes which are connected only at g 0. 

6.4.2 The proof 

The excitations corresponding to the Bethe ansatz (6.52), above the vacuum |Q) = 
\ZZZZ . . .ZZZZ) are symmetric under su{2\2) ® s-u(2|2), each of the branches which 
starts at the node 4 corresponding to one of the s-u(2|2) components. We can therefore 
restrain ourselves to the study of one component, which corresponds to a s-u(3|2) spin 
chain with states Z, X, Y and ip^, ip"^. This spin chain is dynamic, that is its length can 
vary, since the Hamiltonian connects states containing the following two sequences with 
different lenghts [26] 

\...XYZ...) ^ ■ (6.55) 

Terms of this type appear in the Hamiltonian first at order g^, which corresponds to loop 
order one and half at they do not exist at one-loop order. 

The definition of a magnon is, at least at one loop order, a replacement of Z by 
another field X, Y, -0^ or -0^, which are interchanged by the action of su{2\2). At 
higher loop order, because of the dynamical nature of the spin chain, the definition of 
a magnon is more subtle. It is rather easy to define the action of a su{2\l) algebra on 
the magnons, which respects the structure of the interaction (6.55) and does not lead to 
states of fluctuating length. One considers the states \4>) = \X) and |x) = \YZ) together 
with the fermions which form a four-dimensional representation representation of 

su{2\l). The finite-dimensional representations of sit(2|l) were classified in [180], and the 
four dimensional representation generically depend on a continuous parameter b which 
is the charge for the m(1) part of the bosonic subalgebra. In such a representation, the 
generators act like [61] 

R = ^I0)(0|-^|X)(X| (6.56) 

Li = s,^\r){r\-ls'^pm{r\ (6.57) 
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and the odd generators act like 



Q' = a\ij'){<P\ + b\x){ij'\ Q' = a\ij'){<p\ - b\x){ij'\ (6.58) 

-^i = c\i^'){x\ + d\<P){i^'\ S2 = -c\i^'){x\ + d\<P){i^'\ (6.59) 

Prom the commutation relation 

{Q",Sp} = Li + 6'^^{R + C) (6.60) 
with C a central element^ related to the energy 

C = ^ + 6c = b (6.61) 

we obtain the constraint 

ad-bc^l. (6.62) 

Representations with b of opposite sign are conjugate. The quadratic Casimir of su{2\l) 
takes the value b^ — 1/4 in the representation mentioned above, and it vanishes for b = 
±1/2. At this vahic, the four dimensional representation splits into a three dimensional 
fundamental representation and the trivial representation. As we will see, this happens 
at gf = 0, where the interaction term (6.55) does not exist. 

This representation can be hfted to a four dimensional representation of centrally 
extended su{2\2). The commutation relations of su{2\2) are given by [53] 

[Rt,r] = str-ls^r, (6.63) 

[L°,J,] = 6;r-^Sp\ (6.64) 
{Q:,S'^} = 5lLi + 5'^^Rl + 5l6-^C (6.65) 

where the indices a, b, c and a, (3, 7 take the values 1, 2 and J"^ and J'^ are any operators 
carrying the quantum numbers of one of the two su{2) subalgebras. This super algebra 
admits a central extension via 

{Q:,Q^,}^e^^eabP (6.66) 
{5«,5j} = £„;3£«''X (6.67) 

To define a consistent action of the centrally extended s-u(2|2) on the four states of the 
representation previously considered, we have to work with states with variable length. 
We consider a magnon in a periodic chain of large but otherwise arbitrary length 

— ^ n-l,n,n+l, 

|T) = ^e*"f|....ZZ Z T Z ...) with (6.68) 

n 



^Here it was chosen to separate the generator R from C because C is related to the dilatation operator; 
in fact the commutation relations of su(2|l) can be expressed only in terms of the u{l) generator R + C. 
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Figure 5: Structure of the representation (6.71) of the centrally extended sm(2|2) algebra. 
The links in red (green) represent action of two su{2\l) subalgebras. The dotted links 
disappear at g = 0, leaving sm(2|2) four dimensional representations at constant length. 



Upon adding a new site, the one-magnon state transforms as 

\TZ) ~ |T) \ZT) ~ e-'P\T) (6.69) 
The bosonic subalgebra should not modify the length 

Rt = (6.70) 

while for the action of the fermionic operators there are several possible gradings, the 
most symmetric being [181] 

Q:\^') = abWrz'l^) (6.71) 

P\T) = ab\TZ) 
K\T) = cd\TZ-^) 

In this notation, the su(2 — 1) representation displayed in Fig. 5 is realized by \(p'^Z), 
I'i/'^'^Z^/^), at fixed length. The extra sm(2|2) generators connect copies of this represen- 
tations at different lengths. The grading (6.71) makes manifest another su(2|l), with 
representation spanned by \(I)^Z), \tp^''^Z^^'^). 
The closure of the commutation relation imposes 

ad-bc=l (6.72) 

and the energy is given again by (6.60). 

For states with several magnons, one has to take the tensor product of the one- 
magnon representations above, and impose at the end that 1) the resulting state is 
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invariant by (discrete) translation and 2) the global central charges associated with the 
central extension P and K vanish. These conditions are equivalent if one sets 



auh = a{e'P^ - 1) , 0^4 = /3(e-'^'= - 1) (6.73) 
with the constants a and ^ the same for all the magnons. This relation implies 



afe6feCfe4 = K-^ = 4Q;^sin^y or Ck^ ±-^l + lQg''sm'^-^ (6.74) 

where in the last formula we identified by physical reasons a (5 = g"^. The result is the 
BDS dispersion relation, which arises here uniquely from supergroup analysis. Another 
peculiar feature is that, unlike for the usual spin chains, the momentum of the magnon, 
which is related to the spectral parameter, is now a parameter of the representation. This 
is specific to the super-spin chains, since representation of super-algebras often involve 
continuous charges. 

The image presented above is valid for large values of L. For L finite, P and K cannot 
be central elements anymore, since applying K a finite number of times will reduce the 
length of the chain to zero. Since the centrally extended su{2\2) algebra is a contraction 
of the exceptional superalgebra d{2, 1; e) for e — )■ [53], it is tempting to associate e with 
1/L, although this possibility was not explored in the literature. The structure of the 
magnon representation at finite L is still an open problem. 

It is possible to set [53] 

ak^lk, bk^ ^(a^fc-a;^), = ^^-^ , 4 = -— /3(a;^ - ) , (6.75) 

Ikxl axy. 7fe 



which allows to retrieve 



e^^^ = 4, CkAi^ + K-%] ■ (6-76) 



•^k \ -^k -^k 

The closure relation akdk — bkCk = 1 corresponds to the constraint 

+ , 9' _ g' 



+ —r - - — ^ i . (6.77) 



^k ^k 



The dispersion relation in (6.74) appears with both positive and negative sign, and 
the two signs correspond to conjugate representations. On the parametrization (6.75), 
conjugation corresponds to x^ — )■ g^/x^ and it reverses the sign of the momentum as 
well. This is a sign that the conjugate representations should correspond to antiparticles. 
This observation is essential to implement the crossing symmetry which will determine 
the dressing phase. 

The negative branch cannot be seen in perturbative gauge theory. There is no notion 
of antiparticle in the spin chains, and the correspondence with the Hubbard model in 
section 6.2 is able only to produce the positive branch of the energy. This is related to 
the difficulty of obtaining a spin chain description beyond three loops, where the effects 
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of the dressing phase and therefore of the existence of the antiparticles will manifest 
themselves. In the strong coupling or the string limit, which is described, at least in 
some sectors, by a relativistic field theory, the antiparticles are unavoidable. That is the 
case for example for the su{2) principal chiral model described in section 5.2 where the 
inversion x ^ 1/x interchanges the two symmetries su{2)l and su{2)r. It is interesting 
that the two su{2) sectors cannot be separated in the string limit, as noted by Minahan 
[182], as they are in the perturbative limit. 

The next step is to analyze the scattering of magnons. The two-magnon scattering 
matrix can be defined along the lines in section 2.5 for the Heisenberg spin chain, as the 
ratio of the components of the two-magnon wave function for xi <C X2 and X2 ^ xi, 
where xi and X2 are the positions of the magnons. Unlike in section 2.5, here the magnons 
have several flavors corresponding to one of the two centrally extended su{2\2) algebras. 
The two magnon states live in the tensor product of two representations (6.71) and the 
condition that the two-magnon scattering matrix S12 is invariant under the centrally 
extended su{2\2) algebra means that 

[J^ + 4,Si2]^0 (6.78) 

where represent any of the generators in the equations (6.70) and (6.71) in the 
representations associated to the magnon 1 and 2 and characterized by the parameters 
x^, xf. Beisert [53] solved the symmetry constraint (6.78) for the S matrix, which 
depends on the parameters the spectral parameters. The resulting S matrix is 

uniquely defined, up to a scalar phase factor, and, surprisingly, is related to the Shastry's 
R matrix [168] for the Hubbard model [55] [183]. Although this connection to the 
Hubbard model is seemingly unrelated to that described in section 6.2, the identification 
of the coupling constants is the same. One of the implications of this connection is 
that the symmetry which is behind the integrability of the Hubbard model is s-u(2|2), 
although in the Hubbard model part of this symmetry is broken, presumably by boundary 
conditions. 

The centrally extended sm(2|2) symmetric S-matrix obeys Yang-Baxter equation [53] 
[55] and therefore is compatible with factorized scattering. Diagonalizing it [53] through 
the nested Bethe ansatz procedure reproduces the Beisert-Staudacher equations (6.52). 

The link of the S-matrix with the Yangian was investigated in [184] [185] [186]. The 
Hubbard model itself possesses a Yangian symmetry, discovered by Uglov and Korepin 
[187]. 

6.5 The dressing phase at all loops 

The fact that the S matrix is fixed by symmetry up to a scalar factor is similar to 
the situation encountered in integrable two-dimensional relativistic field theories, where 
the scattering matrix of the physical excitation can be obtained postulating a set of 
properties: 

• symmetry 

[Jl^ + J2^<Si2]=0, 
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• factorized scattering (Yang-Baxter equation) 

5i2»5i3»S23 = S23S1SS12 , 

• unitarity 

• crossing symmetry 

S^2iPi,P2)C^SU-p^,p2)C^' = 1 , (6.79) 
where CiS^^Ci^ stands for charge conjugation on the first particle. 

The knowledge of the S matrix allows to construct the Hilbert space as a representation 
of the so-called Zamolodchikov-Faddeev algebra [188] [189], which is the algebra of the 
operators creating the physical excitations. Arutyunov, Frolov and Zamaklar [58] applied 
this idea to the string sigma model in the light cone gauge. The symmetry of the 
excitations is again centrally extended su{2\2) x su{2\2). Counting the number of physical 
excitations, one obtains eight bosons and eight fermions, which in the plane wave limit 
become identical to the ones from equation (5.9). Up to some subtleties related to the 
definition of the action of the fermionic generators, the symmetry in [58] is the same as the 
one used by Beisert to determine the scattering matrix in the gauge theory. Therefore, 
the scattering matrix in the gauge and string theory are the same, up to a twist. The 
condition of crossing symmetry is the same in the two theories. 

6.5.1 The quantization of the su{2) principal chiral model and the AFS dress- 
ing factor 

Let us consider the simple case of the su{2) principal chiral model, where the scattering 
matrix is given by [190] 

s^2\o) = So{ef-±^ with ^o(^) = ^F(iV^T(^ (6.80) 

and 

S12 = 3^2^12 (6.81) 

where, again, the indices R,L refer to the two copies of su{2). The crossing relation 
reads 

So{e + i/2)So{9 - t/2) = . (6.82) 
The excitations have relativistic dispersion relation, with 

E — mcosh(7r^) , p — msinh(7r^) . (6.83) 
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Diagonalizing the scattering matrix via the nested Bethe ansatz, one obtains 



(6.84) 




where L is the circumference of the system, a, ^ — 1, . . . , L, j,l — 1, . . . , and k,n — 
1, . . . , Jy. In the absence of the magnons of type v, these equations are similar to the 
equations for the Hubbard model (6.32) and (6.33), except for the value of So{6) 

This is the quantum solution of the model presented in section 5.2 and it is an 
interesting question how to obtain the classical hmit of the section 5.2 starting from the 
equations (6.84). This problem was addressed in reference [191]. The classical limit can 
be obtained when m ~ e~'^'^^ — )■ and ~ L — )■ oo. One rescales the rapidities as 
6 = 2gC, and reduces the Bethe equations to the equilibrium equations for a system of 
2D Coulomb charges subject to a potential V{^) 



When g ^ oo, the potential V{^) becomes a box with infinite walls situated at ^ = ±1, 
such that the L rapidities 6a/2g = will be confined to this interval. The algebraic 
curve which can be obtained from the equations (6.84) in the continuum limit has four 
sheets, unlike the curve in section 5.2. However, by putting 



the four sheets of the algebraic curve in ^ reduce to two sheets in the variable x, while 
the cut C, € [—1, 1] transforms into the unit circle x — 1. The inversion map x — >■ 1/x 
maps the u variables into v variables, as expected. 

As already mentioned, the equations (6.84) are similar to the equations for the Hub- 
bard model, where 9a ~ 2gsmqa, except for the form So{9). It is therefore interesting 
to check what happens with the equations if a small number of magnons (of type u, for 
example), are kept and if we eliminate the L particles with rapidity 9a. This was done in 
[192], for the leading order in 1/g. The result obtained was the BDS-like equation (6.48). 
Since the result obtained from the principal chiral model cannot be trusted beyond the 
leading order in 1/g, only the AFS dressing factor (6.50) can be obtained. A similar 
result was obtained by Polchinski and Mann [193] starting from an Osp{2m + 2\2m) 
sigma model. 

6.5.2 Higher loop corrections to the dressing phase 

The first one loop corrections to the AFS phase were obtained by Beisert and Tseytlin 
[48] by a direct computation in the string sigma model, and the complete one loop answer 



y(0 = >Cmcosh(27r^O 



(6.85) 




(6.86) 
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was obtained by Hernandez and Lopez [60]. They extended the AFS ansatz to 



OD OO 



a'^{xk,xi) ^ exp2i^ ^ ^Cr,s{g)(qr{xk)qs{xi) - qs{xk)qr{xi)^ , (6.87) 

r=2 s=r+l 

where the first two terms in the expansion of Cr,s{g) at strong couphng are are 

if r + s is odd and zero otherwise. The result was confirmed by Freyhult and Kristjansen 
[194] and Gromov and Vieira [195] [196]. A systematic procedure for the semiclassical 
quantization of the algebraic curves associated with the string sigma model was given in 
[195] [196]. It is difficult to extend this procedure beyond one loop. 

Therefore, the only available option to obtain the full expression of the dressing 
factor is to implement the crossing symmetry. The main difficulty here is that the 
scattering matrix depends on the rapidities of the two particle separately and not only 
on the rapidity difference. This, in turn, can be related to the dispersion relation of the 
magnons, which is not exactly of the relativistic type 



e(p) = y 1 + 16c/2 sin^ | (6.89) 

It is interesting to know that the dispersion relation (6.89) corresponds to a quantum 
deformation of the Poincare group in two dimensions [179] which is in turn related to 
the existence of the elliptic parametrization of the rapidity (6.44). 

While the relativistic dispersion relation can be parametrized in terms of hyperbolic 
functions 

e{u) — mcoshnu , p{u) — m sink nu , — + (6.90) 

and the particle-antiparticle transformation corresponds to u ^ u + i, the dispersion 
relation (6.89) can be uniformized by 

/V 1 , . pis) , sns s 

e(s) = - — , 4osin^-^ = k- — . (6.91 
^ ^ dns ^ 2 dns ^ ^ 

The change of sign of both the energy and the momentum occur for the translation 
s — > s + 2iK', under which we also have 

x^(s + 2iK') = . (6.92) 

x^{s) 

It is this transformation that was used by Janik [61] to implement the crossing sym- 
metry. He worked with the s-u(2|l) representation of the scattering matrix, where the 
problems of length changing are absent. On the labels of the sm(2|1) representation, this 
transformation changes the representation labelled by (1/2, b) to its conjugate labelled 
by (1/2, — b). In order to implement the crossing relation (6.79), one has to specify the 
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action of the conjugation matrix C which relates the two conjugated representations. 
Once this action is specified, the next step is to solve the equation 

Ci5^*^ (si + 2zK', S2)Cr' = S{^{sr, S2) , (6.93) 

in the tensor product of representations (1/2, bi) and (1/2, b2) . The symbol represents 
supertransposition in the first space and we used the uniformization parameter s for 
the rapidity. For generic values of bi and b2, the tensor product decomposes in three 
irreducible representations, and equating the coefficients of the projectors on these three 
representations allows to obtain in particular the condition satisfied by the scalar part 
of the scattering matrix 



+2«',..)s.(.„..) ^ ggigg::;; . (6.94) 



Similarly to the su{2) case, we do not expect that the dressing factor is invariant by 
shifts of 4iK' in Si. 

The crossing equation (6.94) has many solutions. Among the required properties, 
it should reproduce the known AFS [51] and Hernandez Lopez [60] behavior at strong 
coupling and the behavior at weak coupling 5*0 = 1 + 0{g^) at weak coupling. 

A solution which satisfies these constraints was constructed by Beisert, Hernandez 
and Lopez [62]. This solution can be written in the form (6.87), 

S-^{xuX2) ^ I i: V/T- ^"'(^i^^2), (6.95) 

oo oo 

a{xi,X2) = expz^ Y g''+^-'^Cr,s{g)(qr{xi)qs{x2) - qs{xi)qr{x2)^ , 

r=2 s=r+l 

with the coefficients Cr^s{g), non zero only for r + s odd, expressed as a series in 1/g 

oo 

c.,3(^) = ^4:)<7-" (6.96) 

n=0 

and the numbers ci^] are determined as 

(„) _ C(n)(r-l)(.-l)r[l(r + . + n-3)]r[i(.-r + n-l)] 
""^'^ (-27r)-r(n - 1) r[l{r + s-n + 3)]r[lis-r-n + l)r^'^''^'"'' ^''•^'^ 

For n = 0, 1 the coefficients Cr"i are defined as in (6.87). Beisert, Eden and Staudacher 
[63] gave compelling evidence that the coefficients describing the dressing function can 
be expressed by an expansion for small g in the remarkable form 



oo 



^(^) = E"^7^^" (6.98) 



n=0 
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where Cr,s" can be defined for positive n by analytical continuation of the ones in (6.97) 

= (6.99) 

C(l + n) cos [f ] (r - 1)(5 - l)r(n + 2)r(n + 1) 6r+s,i mod 2 

r [i(5 + n - r - s)] r [i(3 + n + r - s)] r [i(3 + n - r + s)] r [i(l + n + r + s)] ■ 

As expected, the pertubative expansion contains only even powers of g, due to the factor 
cos7rn/2. As ({n) — > 1 for n — > 00, (6.98) is convergent around g = 0, while the 
series (6.96) in 1/g has zero radius of convergence. This fact has deep consequences and 
enormously complicates the strong coupling limit. The analysis of the strong coupling 
limit is the subject of the next section. 

Being able to express the dressing factor in perturbation in g allowed to check that 
it is trivial up to three loop order and to compute anomalous dimensions recursively up 
to the desired order, at least for operators without wrapping corrections. In particular, 
Beisert, Eden and Staudacher [63] computed the four loop contribution to the so-called 
cusp anomalous dimension. Remarkably, this correction coincided with an impressive 
four loop computation [64] of the cusp anomalous dimension from the gluon amplitudes, 
later confirmed in [197]. The four-loop direct computation of the dilatation operator in 
the su{2) sector [141] confirmed the result. 



6.5.3 The "magic" BES formula 

Since the dressing factor is an important ingredient for determining the anomalous di- 
mensions, it is highly desirable to express it in a more compact form than that in (6.96) 
and (6.97). A remarkable representation as a convolution of scattering kernels was given 
in [63]. Let us define the scattering kernel, as in appendix B, as the derivative of the 
scattering phase 

K{u,u')^^^ipiu,u') (6.100) 
and take its Fourier transform on both variables u and with the definition 



/OO POO 
du / du' K{u,u') e^(*"+*'"') . 
-00 J —00 

An important building block of the Bethe equations is the sti(l|l) kernel 

Km{u,u') = -2X5„(i|i)(li,u') = 



(6.101) 



(6.102) 



i d 
TT du 



In 1 - 



- In 1 - 



X (u)x'^(u') 



x'^(u)x (u') 

After Fourier transform, due to the structure of the singularities in u and u', we obtain^ 
K^{t,t') = 8ng\l - signtt')KI ^^:M^:M^e-(l*l+l*'l)/^ . (6.103) 



n>0 



^Our definitions for tlie kernels differs from tlie ones of [63]. We try to keep the notation unitary and 
as close as possible to the different notations in the literature. 
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To avoid complications due to the sign of t, we define tfie Fourier transform only for 
positive values oi t,t'. In the next chapter, we are going to explain the significance of 
this simphfication. We define then 

K^{t,t') = 167rgHJ2n ^ff^ff e-('^''y' ^ K+{t,t') + K.^ , (6.104) 

n>0 ^ ^ 

where and K_ contain the expansion on odd and even Bessel functions respectively. 
The "magic" formula of BES expresses the dressing kernel as 

K4t,t') = —K.{tX)Y:r^K4t",t') (6.105) 

where Kd{t, t') is defined as 

• /"OO poo 1 

Kd{t,t')^- du cia'-f ln(7(ii,tiO e^(*"+*'"') . (6.106) 

TT J-oo du 

Inverse Fourier transformed versions of the dressing phase appeared in [78] and [65]. The 
Dorey, Hofman and Maldacena formula 

d{x,y) = x{x'^ ,y^) + x{x~ ,y') - x{x^ ,y~) - x{x~ ,y^) (6.107) 

/ \ . f dzi f dz2 1 , T^/. • / 1 1 \\ 

X{x,y) = f TT f ^7 V7 r lnT{l + tg{zi + Z2 ). 

J 2n J 2-K {x - zi){y - Z2) zi z^ 

was used to prove [66] [67] that the dressing factor satisfies the crossing equation (6.94) 
for arbitary values of the coupling constant g. 

An attempt to derive the formula from the consistency conditions of the ps-y,(2,2|4) 
Bethe equations was made by Sakai and Satoh [198] [199], although the interpretation 
of their calculation is still mysterious. Before, a formula with a similar structure was 
obtained for some special state inside psu{2, 2|4) [200]. Janik and Lukowski [201] argued 
that this pattern is rather common among the 2d integrable field theories. 

7 Strong coupling limit of the asymptotic Bethe An- 
satz equations 

7.1 The cusp anomalous dimension 

The most sophisticated checks of the Bethe ansatz equations (6.52) involve the anomalous 
dimension for the twist two operator, or the so-called cusp anomalous dimension. In the 
language of the J\f — A SYM theory, the twist two operator can be defined as 

Tr{ZV^Z) + ... (7.1) 

where I) = + iVi is a light-cone component of the covariant derivative, the dots 
represent terms with different orderings of the fields, and S is the Lorentz spin. Gross 
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and Wilczek [202] showed that the corresponding QCD operators has, at one loop in 
perturbation theory, a dimension which scales logarithmically with S. Later [203] [204] , 
arguments were given that the logarithmic scaling should hold at any order in pertur- 
bation theory. Korchemsky [203] showed that the anomalous dimension of the twist-two 
operator is twice the cusp anomalous dimension, which appears in the divergences asso- 
ciated to the Wilson loops with cusps [68]. Therefore, the conformal dimension of the 
operator (7.1) scales at large S as 

A = S + 2 + f{g)\nS + ... (7.2) 

where f{g) will be referred to in the following as the universal scaling function. The 
denomination of universal originates in the fact that this factor does not depend on the 
value of the twist, at least for finite value of the twist. A twist L operator is defined as 

Ti {V'^'ZV'"' Z ...V''^ Z) + ... with h + k2 + ...kL = S (7.3) 

The logarithmic scaling persists when the twist scales like the log of the Lorentz spin, 
L ~ In^ [205], although the scaling function depends now on the new parameter, see 
also [87] 

J = j^, f{9)^f{9:j)- (7.4) 

The universal scaling function was determined in perturbation in the gauge theory, at 
three loop order in [206], who extracted it from a computations for QCD [207], and to 
the fourth order numerically from the four-gluon amplitude in [64], [197], 

m = 8/ - + gvrV - 16 (^vr^ + 4C(3)^) / + 0(^^°) (7.5) 

where we have replaced the numerical value [64], [197] by the probable exact value. 
From the string point of view, the state (7.1) corresponds to a solution for a folded 
string spinning in AdS^ [72]. The leading term in the strong coupling expansion was 
computed in [72] , while the one-loop order was computed in [73] and the two loop order 
in [74, 75] 

where K is the Catalan constant. The two expansions (7.5) and (7.6) were reproduced 
from the Bethe ansatz equation and they can be performed recursively to arbitrary order 
in g or 1/ g. This result can be obtained by solving an integral equation for the magnon 
density corresponding to the state (7.1) with lowest energy, which belongs to the rank 
one sector sl{2), the non-compact version of su{2). The resulting integral equation is 

p(«) +P/j(«) = -2:^^ + y du' )C{u,u') p{u') , (7.7) 

where p{u) and Ph{u) are the densities of magnons and holes respectively, normalized 
such that 

J du p{u) = S , J du Ph{u) ~ L (7.8) 
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and the kernel lC{u, u') is composed from three parts 



K,{u, u') = KsH2) {u, u') - K^{u, u') - Kd{u, u') . (7.9) 

The kernels Km{u,u') and Kd{u,u') were defined in (6.102) and (6.105) and they vanish 
for = 0, while 

Ksl(2) («, u') = -Ksu{2) («, u') = ^ _ ^ ^ . (7. 10) 

In the lowest energy state, L — 2 holes are closely packed around u — and two are 
situated at around u = S/\/2, see e.g. [84, 208]. The equation (7.7) was solved, for L/S 
finite and g = 0, in [209] where it was shown that the density p{u) with u = u/S has 
supports on two intervals on the real axis symmetric with respect to the origin, [—6, —a] 
and [a, 6]. When L/S 0, the gap around the origin closes, a — )■ 0. The way the gap 
closes depends on the way L/S goes to [205]. When L is finite and g — the equation 
(7.7) becomes singular. Korchemsky [210] solved this equation in the rescaled variable 
u — u/S 

1 l + y/T^^ 
Pq[u) = - In , . 7.11 

TT 1 - V 1 - 4^2 

Let us note that the density is normalized as J du poiu) = 1, which expresses the 
fact that the total umber of magnons is 5". At S — > the Korchemsky density has a 
logarithmic divergency ^0(^0 ^ ~^ ln(2'u^). This divergency can be removed if one solves 
the equation in the original variable u around u — [210, 211, 208] to get 



Po{u) ~ -- 

TT 



2 

V^(l/2 - iu) + V^(l/2 + iu) + In 



(7.12) 



At scales 1 < u ^ S, the leading contribution in^ in po{u) is the constant part po = 
I In S. Eden and Staudacher [70] proposed to treat the case 51 7^ by separating the one- 
loop contribution from the rest and writing an equation for the density for the fiuctuation 
a{u) defined by 

p{u) = [ao - a{u)] \nS + 0{S^) (7.13) 

with (To = ^. This separation works provided we take the limit 5 — )■ 00 and neglect 1/S 
corrections. Then a{u) — > when u goes to infinity. The equation satisfied by a{u) is 
called the BES equation [63] 

/oo 
du' ]C{u, u') {a{u') - (7o) (7.14) 
-00 

and the universal scaling function is, up to a constant, the normalization of the density 

a{u) 

POO 

f{g) = 2 du a{u) (7.15) 
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This expression, due to Lipatov, can be obtained using the expression for the energy 
and the integral equation (7.14). The fact that a{u) has a non-zero norm may seem 
paradoxical, since it is the leading term in the expansion in 1/S of difference of two 
densities po{u) and p{u) which are normalized to S. The reason is that the expansion 
in l/S* of the density is not uniform in u. Another remark is that a shift of the density 
a{u) by a constant leaves the equation (7.14) unchanged. Therefore, the right solution 
is obtained from the condition a{u) ^ when u ^ oo, which is necessary in order to 
have a finite normalization. 

The BES equation can be solved relatively easily in perturbation in g hj a recursive 
procedure, since the expansion of the kernel K{u, u') around g is well-behaved. The first 
four orders in g"^ obtained in [63] coincide with the ones obtained in perturbative gauge 
theory (7.5). The expansion m. 1/g was much more difficult to obtain, due to the fact 
that 1/(7 = is an accumulation point of singularities. The limit of the kernel K(u,u') 
is not uniform in u and one has to deal differently with different regions in the rapidity 
space u [78, 212] . The leading term in (7.6) was first obtained numerically [71], and then 
analytically in [76, 77, 78, 79]. The next term was obtained by Casteill and Kristjansen 
[80] and Belitsky [81], and the third term and a recursive procedure for the next terms 
by Basso, Korchcmsky and Kotanski [82] and in [83]. 

The finite-twist operators in A/" = 4 SYM were extensively studied, including the 
wrapping corrections at five loops for the twist-three operator [213] and logarithmic 
corrections to the scaling with the Lorenz spin [214]. The finite twist operators obey 
an interesting property called reciprocity [215] [216], which was observed previously in 
QCD [217]. Reciprocity was studied from the point of view of Bethe ansatz in [218]. 

7.2 The generalized scaling function 

Essentially the same methods which were used to compute the universal scaling function 
can be used to analyze the generalized scaling function f{g,j). Freyhult, Rej and Stau- 
dacher [84] proposed use parametrization similar to (7.13) but now for the ensemble of 
magnon and hole densities 



The support of the hole density is the interval [—a, a] and the support of the magnon 
density is the complement M\[— a, a]. The density fiuctuation is now subject to the FRS 
equation [84] 



p{u) + ph{u) - [ao - a{u)] ln5 + 0(5°) 



(7.16) 




(7.17) 



The relation between j and a can be obtained via the normalization condition 




(7.18) 



and the generalized scaling function is given by 




(7.19) 
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At small g and j, the equation (7.14) can be solved perturbatively in both parameters. 
The equation (7.14) was rederived and analyzed perturbatively in j at any g up to eight 
order in j in [85]. When a <^ 2g, the hole density can be approximated with the density 
of excitations in the 0(6) sigma model [86]. In this regime 

jr^m = kg'/'e--%l + 0{l/g)), with ^ = (7-20) 

The relation to the 0(6) sigma model was previously identified in [87]. A different regime 
corresponds to a > 2g, where the relevant parameter is 

(7.21) 

In this regime, the universal scaling function was shown to behave [88] at strong coupling, 
for small £ and up to two loop order 

f{9,j) = f{9)-J (7.22) 



+ f 



2, + 1 Q - 21n£) + ^ (go2 - 61n£ + 81n^£) + (1) 



+ o(r) 



which means that the i = point is singular, therefore the limits j — and (7 — )■ oo do 
not commute. This structure can be predicted from the connection to the 0(6) sigma 
model [87]. The coefficient go2 was predicted from an early string computation [88] to be 

g?r^ = -2K-^ln2 + ^ (7.23) 

while the computation [89] using the Bethe Ansatz equations, confirmed by [90, 91] 
predicts 

5o^,^ = -^ln2 + ^. (7.24) 

Recently, a revised string computation [92] reproduced the Bethe ansatz result (7.24). 

Finally, it would be interesting to analyze the transition between the two previous 
regimes, especially the regime when a — 2g 1, which should occur for values j ~ g"^^^- 

Other computations in the strong couphng limit concern the correction [219] [220] to 
the logarithmic scahng, the so-called virtuality BL{g), defined as 

A - S - L = f{g){\nS + 1e-{L- 2) ln2) + BL{g) + 0{l/S) (7.25) 

where 7^ is the Euler-Mascheroni constant. 



7.3 The strong coupling limit 

Several approaches were taken to approach the strong coupling limit of the BES equation 
(7.14) or other similar integral equations. One of them proposed in [71] and which was 
successfully used for numerical computations, uses the representation of the Fourier 



66 



transforms of the kernels in terms of Bessel functions (6.104), which are used as building 
blocks of the magnon densities. This method was used in [71, 77, 82]. Other authors used 
the rapidity representation, starting either from the BES equation, like in [76, 78, 83], 
or directly from the Beisert-Staudacher equations for large twist L and then taking the 
small L limit, like in [80, 81]. 

Here, the method of [83] will be presented. It allows to obtain the original result of 
[82] by using standard analyticity assumptions with a transparent physical interpretation, 
and in the same time it may be helpful to make contact with other problems at strong 
coupling, like to quantization of the algebraic curve (4.29). 

In order to avoid clumsy formulas, we are going to use notations which are better 
adapted for the strong coupling limit 

— '"old /„ r,a\ 

e=-, — , (7.26) 

as well as the variable x{u), related to u by 



Note the branch cut of x{u) for m e [— 1, 1]. In the intermediate steps we will also use 
the notations 

x^{u) = x{u±ie). (7.28) 
The kernels will be rescaled by Ag^ and the densities by 2g. 

7.3.1 The resolvent 

Instead of solving the BES equation for the density, we are going to write the equation 
satisfied by the resolvent 

i?phys(w)= / du' (7.29) 

The resolvent is analytic in the complex plane except on the support of the density. If 
the density is normalizable, the resolvent behaves at — > oo as 

i?phys(M) - - / du' a{u'). (7.30) 
In the case we are going to consider, the support of the density is the real axis, so that 

= ;r- [Rphysiu - iO) - -Rphys(^i + iO)] (7.31) 
The resolvent is closely related to the Fourier transform of the density, cr(t), defined as 

/•oo 

a{t) = du e^*"a(M) . (7.32) 
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Assuming that u is in the upper half plane (UHP), we have 

/oo poo 

du' / dt e^*("-"'V(ii') (7.33) 
oo Jq 
/•oo 

= -i / dt e^'^'ai-t) . (7.34) 



Since the density we are working with is symmetric a{u) — a{—u), the resolvent is 
antisymmetric i?phys('i*) = ~-Rphys(~'"); so we can retrieve its value in the lower half plane 
(LHP) once we know that in the UHP. In the following we are going to call resolvent the 
function R{u) = i?phys(^) in the UHP and which is analytically continued to the LHP. 
Therefore, R{u) is analytic in the UHP, but it may have singularities in the LHP and 

a{u) = — — [R{u + iO) + R{-u + iO)] (7.35) 
7.3.2 The linearized BES equation 

Our starting point is the BES equation written in Fourier transformed form 

df' 

[1 - i^.U2)W] ^(^) = - ^ {K^{t,t') + Ka{t,t')) {a{t') - ao(t')) , (7.36) 

where ao(t) = Se5{t) and Ksi(2)(t) — e~^^*. In abstract notations it can be written as 

-5" V = {K+ + K_ + 2K_SK+) {a - (Tq) , (7.37) 

where we used 

S-' = l- K,i^2) , (7.38) 

with Kj^ and K_ defined in (6.104) and the dressing kernel is given by the "magic 
formula" (6.105). The relation (7.37) can be rewritten as 

-2(7 = [(1 + 2SK+){1 + 2SK_) - 1] (a - ao) , (7.39) 

This equation can be transformed into a system of equations where the kernels K± appear 
linearly, at the expense of introducing an extra density r such that 

T + a = -2SK+{a-ao) (7.40) 
T-a = -2SK_{t + ao) . 

This procedure was used in different forms by Kotikov and Lipatov [76], Eden [221] and 
Basso, Korchemsky and Kotanski [82] and is essential to obtain equations which can be 
exploited at strong coupling. Going back to the rapidity representation, these equations 
can be written as 

H{u + 2ie) - H{u) + R{u + 2ie) - R{u) = 8ie- — - + 2 du' K+{u, u')R{u') 

dux'^ J 

H{u + 2ie) - H{u) - R{u + 2ie) + R{u) = 2 J du' K^{u, u')H{u') (7.41) 

68 



where H{u) is an auxiliary resolvent associated to the density t{u). These equations are 
derived for > and then analytically continued for any u. A set of similar, equivalent 
equations hold for < 0. Since R{u) and H{u) are analytic in the UHP, any term in 
K±[u, u') which is analytic in u' in the UHP will give a zero contribution to the integrals 
in the r.h.s. We therefore can set without loss of generality 
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V ^' y 




1 d 


In 


\- 1 




27ii du 


\ x+y- 



+ In 1 + 



In 1 + 



x^y^ 



(7.42) 



where we used the notation 
u±ie = - 



x^ + 



X 



-t 



u' ± ie 



y 



The key observation is that the kernels K±{u,u') depend on the coupling constant e = 
l/4gf only through the shifts u,u' ± ie, and the shifts can be absorbed by a change of 
variable and the shift of the integration contour. The equations can be rewritten in terms 
of the kernels 



lim Km(u, u') , 

e-i-O ^ ' 



(7.43) 



27ri 1 — x"^ \y — - v -\- - 



which depend on the variables x, y with a cut on [—1, 1] with x{u + iO) = l/x{u — iO). 
Upon introducing 



R±{u)^-[R{u)±H{u)] 

as well as the related functions r±{u) 

r±{u) — R±{u — ie) — R±{u + ie) . 
we arrive at the equations 

r+{u) = r - / du' K+{u,u')[R+{u' + ie) + R^{u' + ie)] 

X — 1 JR-iO 

-- / du' K_ {u, u') [R+ {u' + ie) - R_ {u' + ie)] . 

jR-iO 



(7.44) 



(7.45) 



(7.46) 



The r.h.s. of these equations have a single branch cut of order two for m e [—1, 1] and they 
have definite symmetry properties under u — > —u. The same properties are shared by 
the l.h.s., so we deduce that r^{u) and r^{u) are symmetric, respectively antisymmetric 
functions with a single branch cut on [—1, 1]. At infinity they behave as 



r_|_(M) oc 



u^ 



r-{u) oc 



{u oo) 



(7.47) 
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and at — > ±1 they have singularities at most of the type 



The relation (7.45) can be inverted to obtain 

CO 

R^{u) = r±{u + {2n + l)ie) . (7.49) 

n=0 

Therefore, the resolvents R±{u) have a set of equidistant branch cuts in the LHP, sepa- 
rated by a distance 2ie. In the strong coupling limit, the distance between cuts vanishes 
and we are left with a single cut on the real axis plus non-perturbative terms. The con- 
tour of integration in (7.46) gets pinched between two branch cuts. It is therefore better 
to define the action of the kernels K± in a manner which can be exploited at strong 
coupling. 

7.3.3 The functional equations 

On functions F{u') analytic in the UHP with the real axis included (as R±{u' + ie) are), 
we can deform the integration contours for u' defining the action of the kernels K±(-u, u') 
such that the new contour encircles the branch cut u' G [—1, 1], and then transform the 
contour integral into a line integral using y{u' -|- iO) = l/y{u' — iO), 

l-x^ '27rt\y-x y + x y-l^y+l) ^> 

= / T -, T -, Fiu') . (7.50 

J_^+i^ 2t:i x-l\u'-u^ u' + u) ^> ^ ' 

Finally, the second term in the last line can be transformed, using u' — > —u' , so that we 
get 



In this form, the kernel can be applied to functions with a branch cut on the real 

the resolvents H{u) and R{u) become in the limit e = l/4g( 0. The kernel in equation 

(7.51) is very similar to a Cauchy kernel. If one defines the function F±{u) by 

\±F[—U) if iSU < , 

then the kernels K± can be written as Cauchy integrals 



K±F(«) = (t^ F±{u') . . (7.53) 

^ ^ ^ 27n ^ ^ \ u^ -1 u' -u ^ ' 
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This definition is useful when F±{u) is analytic outside the [—1, 1], in particular when 
F±{u) — F{u), which means that F{u) is symmetric/antisymmetric. On the space jC± 
of functions f± which arc i) symmetric/antisymmetric, ii) analytic everywhere except a 
branch cut on [—1, 1] and iii) vanish faster than 1/u at infinity the kernels K± act as 

K±/± = /± . (7.54) 

Since K±f{u) belongs to the space C±, for any function f{u) we deduce that the kernels 
K± are idempotents 

K| = K± . (7.55) 

Finally, we notice that for the constant function f{u) = 1, which does not belong to the 
spaces jC±, we have 

1 = , X_ 1 = . (7.56) 

1 — 

At this point, we can express the equations (7.46) as an equation for the zero modes of 
the kernels K±. By redefining 

r+ (u) + r_ (u) = R+(u - ie) + R-(u + ie) + Aie 

T+{u)-T^{u) = R-(u - ie) - R+(u + ie) + Aie . (7.57) 

we arrive at the functional form of the BES equation 

K+(r+ + r_) = o, (7.58) 
K_(r+-r_) = o. (7.59) 

As it can be easily deduced from the expression (7.51), the zero mode condition can be 
expressed as a functional equation on the cut. The two equations (7.58) are equivalent 
to the functional equation 

r+{u + iO) + T_{-u + iO) = 0, we [-1,1]. (7.60) 

Note that this functional equation is valid at any value of g, and it is equivalent with the 
original BES equation, provided we specify the class of functions to which T±{u) belong, 
namely functions with a collection of equidistant square root cuts 

11°° 
r^u) = -r+{u) - -r_{u) + ^r+(M + 2me) , (7.61) 

n=l 

I 1 °° 

T+{u) = Aie + -r+{u) + -r_{u) + ^ r_(M + 2nie) . (7.62) 



n=l 
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7.3.4 The perturbative solution at strong coupling 

At strong coupling, the distance between the square root cuts vanishes and the func- 
tions T±{u) reduce to a functions with a single cut on [—1, 1] and a non-perturbative 
contribution. Let us consider the combinations 

r±(M) Tr±(-'u) = T^r^(M + 2me) (7.63) 

The functions in the r.h.s are periodic in u with period 2ie, therefore they are purely 
non-perturbative. In the perturbation expansion in e they can be put to zero, so that 
the functions ^±{u) are now symmetric/antisymmetric 

r±H = ±r±(-«) . (7.64) 

Introducing the functions 

G±{u) = T+{u) ± iT_{u) . (7.65) 
the symmetry property becomes (7.64) becomes 

G+{u) = G_{~u) . (7.66) 
while the functional equation (7.60) becomes 

G±{u + iO) = ^iG^{-u + iO) = TiG±{u - iO) . (7.67) 

This equation can be considered as a monodromy equation for G±{u) around the cut. 
The monodromy is of order four, so we deduce that G±{u) have a branch cut of order 
four between [—1, 1]. Introducing the uniformisation variable s via 

u-1 ^ ^ 

we can write the general solution of the monodromy equations (7.67) as 

G±{s) = 4ze^c„(e)e±(2"+i/2)^ 
nez 

^ , (7.69) 

where the coefficients c„(e) are to be fixed. Basso, Korchemsky and Kotanski [82] sup- 
posed the coefficients c„(e) can be expanded in a series in integer powers in e. We suppose 
the same thing here, although this is not generally true for generic states, see for example 
the su(l|l) sector analyzed in [222]. Prom the behavior at infinity of the two functions 

n>0 ^ 
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we can extract the universal scaling function, which is encoded into the physical resolvent 
Rphys{u), according to equation (7.30). We obtain the conditions 

The coefficients c„(e) can be determined, at least recursively, by analyzing the struc- 
ture of the general solution (7.69) in the vicinity of the points u — ±1. We introduce 
the rescaled variables 

«-l _ u + 1 
z^^, . = (7.72) 

The variable z coincides, up to a shift by 2g, with the original (before rescaling by 2g) 
rapidity in the BES equations. In the variable z, the functions G± still have a set of 
distinct cuts separated by i. When g ^ oo, these cuts become semi- infinite, extended to 
from 2; = 0tO2; = — oo, and they can be expressed in terms of the functions r±{z) which 
have, as discussed in section (7.3.2), the following expansion 

r±{z) = b^{e) z-'/' + Y dni^) ^" ( kl < 1 )• (7-73) 

n>0 n>0 

For a more symmetric presentation, one can introduce the notations 

g±{z)^r4z)±ir4z), (7.74) 
The functions g±{z) are related to G±{z) via 

g±-^^iD'-l)G^, (7.75) 

with D the shift operator Df{z) = f{z + i/2). 

The compatibility of the expansion (7.73), and the similar expansion around z = 0, 
with (7.69) requires that 

oo 

c„(e) = el"l an{e), an{e) = <^n,p {ne'L). (7.76) 

p=0 

In particular, this implies that at the leading order in e 

00 = 1 + 0(6), and /(^) = ^ + 0(6°) 
with the corresponding resolvent 

which is the solution which was obtained in [77] [78]. 
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In order to determine the higher orders in e, we have to implement the compatibihty 
condition between the (large) u expansion (7.69) and the (small) z expansion (7.73). We 
first write 

and then take the inverse Laplace transform defined by 

f(^) = ^if dze^'fiz). (7.79) 



iR+O 



The explicit expression for the inverse Laplace transform of (7.78) is a series of confluent 
hypergeometric functions of the first kind 



nez \ / V / 



(7.80) 



where in the second line we have used the transformation properties of the hypergeo- 
metric function and we have neglected terms of the order e"^^^. This is an expansion in 
quarter powers of i, while g±{i) will contain half- integer powers of £. Under the inverse 
Laplace transform, the relation (7.75) becomes 

fcM = ^^5.(q. (T.81) 

sm(| ± f ) 

Now we represent, as in [82], the ratio of the sine functions in (7.81) as 

sin(l) _ s±{e) 



sin(| ± f ) T^{e) ' 
where S and T represent ratios of Gamma functions: 



(7.82) 



^\2tt) ^ ^2 27r ^ 4^ 

which expand at ^ = -|-oo and at ^ = respectively as 

S±{i) = ± (£/27r)^^^ 1 + J2s^e-^ ) , (7.84) 



n=l 
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If we rewrite the equation (7.81) as 

G±{i) ^ 1 9±i^) 



(7.86) 



then the l.h.s. is analytic everywhere except the negative real axis, while the r.h.s. 
is analytic everywhere except the positive real axis. As a consequence, neither of the 
sides has poles and the only singularities can be branch points at £ = and £ = oo. 
This means, in particular, that the expansion of the r.h.s. as a power series at ^ = oo 
coincides with the expansion of the l.h.s. at ^ = 0. Both sides expand in the same 
fractional powers, so that if we multiply both sides of equation (7.86) by by (£/e)^^^/^ 
we obtain 

^{£/er'^ =E CHe)r\ (7.87) 
where the coefficients C^(e) should be understood as formal series in e, 

oo 
p=0 

Prom (7.73), (7.84) and the relation (7.86) we deduce that the coefficients in front of the 
non-negative powers of i vanish, 

C^(e) = for n = -l,-2,... . (7.89) 

Solving the contraints (7.71) and (7.89) order by order in e one can evaluate recursively 
the Taylor coefficients of the series (7.76). The recurrence procedure is possible 
because at each order in e the sum on the r.h.s. of (7.87) contains only a finite number 
of negative powers of i. 

In the leading order in e the series expansion of G± at £ = 0, given by 

G±(£) = 2i(6/£)^^i f t + 0{e) \ , (7.90) 

contains only non-negative powers in £. Therefore the sum on the r.h.s. of (7.87) contains 
only the term with n = 0, and we have 

oo 

(£/6)^^iG±(£) = 2. J]-^±^r =r±(£)Co±o- (7.91) 

n=0 M"'=^ 4^ 

Prom the constraint (7.71), which in the leading order gives q;o,o = 1, we evaluate 

^0^0 = 2^^, Co;o = 2i^^. (7.92) 
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For the other coefficients we find 



where are the coefficients in the expansion (7.85). Continuing this procedure order 
by order in e, we reproduce the result of [82] for the universal scaling function, 

/(e) = - + 4j]eH-ina„(e) (7.94) 

n=l 

1 31og2 K 
= ^ ^e + ... . 

The method presented here can be extended to the computation of other quantities, 
in particular it can be adapted to the computation of the generalized scaling function, as 
it was done by Volin in the regime j ^ g [91]. It was also used by Volin [223] to compute 
analytically the mass gap for the 0{N) sigma model. 



8 Finite size corrections and the thermodynamical 
Bethe ansatz 

In the current state of the development of the subject, there is compelling evidence 
that the dilatation operator of the M — ^ SYM theory is integrable and the results for 
its spectrum are compatible with the AdS/CFT correspondence. The above picture is 
coherent, although the problem of computing conformal dimensions for operators of finite 
length is not yet solved. The Beisert-Staudacher equations (6.52) are correct for lengths 
of the chain L large. The finite size corrections are associated to the so-called wrapping 
interactions because they involve spin interactions that wrap around the circumference 
of the spin chain. Those terms appear in perturbation theory as terms of the order g"^^, 
hence they are exponentially suppressed with the size of the chain. A first glimpse of 
how they should appear in a spin chain is given by the Inozemtsev model. There, the 
Bethe ansatz disappear at finite L and the two-magnon energy has a non-additive part 
which also vanishes exponentially with the system size [33]. Prom an (integrable) field 
theory point of view, the finite size corrections shoTild come from the virtual particles 
propagating around the cylinder [93]. A method to compute these terms for a generic 
two dimensional field theory was given by Liischer [95, 96]. The Liischer method was 
implemented for AdS/CFT by Janik and Lukowsky [94] and it was pursued in a vast 
number of circumstances, culminating with the computation of the wrapping correction 
at four [97] and five [98] loops for the Konishi operator . The result was shown to agree 
with the diagrammatics computation in A/" = 4 SYM [99, 100]. The same method was 
used to show [224] that the wrapping corrections are compatible with the analyticity 
requirements imposed by the BFKL equation on the anomalous dimension of the twist 
two operator [225] and which were not met by the asymptotic Bethe ansatz [226] . 

A method to compute the finite-size corrections for integrable field theories and inte- 
grable models in general is the so-called Thermodynamic Bethe Ansatz (TBA). Thermo- 
dynamic functions for integrable systems where first computed by Yang and Yang [227] 
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for the Bose gas with delta interaction and soon after by Gaudin for the XXZ model 
[228] and Takahashi [229] for the XXX model. The distributions of roots were found 
to obey certain functional equations, known in recent years under the name of Y sys- 
tems. Similar relation were found around 1990 among the transfer matrices for various 
integrable models [230, 231, 232]. These relations, also known as fusion relations, can 
be encoded in terms of the Hirota equation. For models based on Lie supergroups of 
type su{M\N) Hirota equation was first obtained by Tsuboi [233]. A connection be- 
tween these seemingly different equations, at the origin of the method for computing 
finite size corrections using the thermodynamic Bethe ansatz [234] , is the interpretation 
of the finite size L of an integrable system as the inverse temperature in the system 
with the space and time interchanged. The system of functional equations associated 
to the thermodynamical Bethe ansatz contains data from a Dynkin diagram [235]. An 
equivalent formulation of TBA in terms of a non-linear integral equation was given by 
Destri and De Vega [236]. The TBA approach was initially used to compute the finite 
size corrections to the vacuum energy, but later it has been understood how to extend 
it to excited states [237, 238] . 

In a relativistically invariant system, the rotated system, called mirror theory in [102], 
is equivalent to the initial one. In the context of AdS/CFT, because of the choice of 
the light cone gauge, the excitations do not have a relativistically invariant dispersion 
relation, but the BDS dispersion relation (6.3). The systematic study of the mirror model, 
including the bound states, was initiated in [102] and continued in [239, 240, 241]. 

The TBA method was recently applied to the AdS/CFT system [104], after working 
out in detail the simpler su{2) principal chiral model [103], with which it bears some 
similarity. In particular, this method allowed to compute the four loop wrapping correc- 
tion to the Konishi operator [104] and the five loop numerically [105] and analytically At 
strong coupling, a numerical analysis of the Y-system was done [242] which confirmed 
the leading A^/^ and permitted to fix the first few coefficients in the strong coupling 
expansion. Compared with the prediction from the string theory [243], the result differs 
in the coefficient of the X^^^^ term. In [244, 245] the finite size corrections for operators 
with L M g were analyzed and the corresponding results were compared with the 
predictions from the quasi-classical analysis of the algebraic curve. 

In [246, 247] the Y-system was derived from the Beisert-Staudacher equations for the 
mirror model. The Y-system is a system of finite-difference equations obeyed by the 
densities of magnons/holes at finite temperature. These equations can be restated as 
Hirota equations [103, 104], which are usually obeyed by the transfer matrices associated 
to particular representations of su{n) [248] or su{n\m) [249] algebras. 

In the following, I will explain how the formalism of the Hirota equation and Y-system 
can be defined for the su{2) principal chiral model, which is a simpler version of the Y- 
system which appears in the AdS/CFT problem. Then, 1 will present the AdS/CFT 
Y-system, which was used to compute the dimension of the Konishi operator at weak 
and strong coupling. The presentation follows mainly the references [103, 104]. 



77 



8.1 Y system for the principal chiral field 

We are starting with the Bethe equations for the principal chiral field in a large volume. 

- 11 s,{e^ - e,) 11 11 ^^_,,_,/2 ' 

with the scattering matrix 5*0 (6*) satisfying the crossing equation 

S,ie + t/2)So{e-z/2) = ^^^. (8.2) 
The energy of a state is given by 

E = m cosh TiOg . (8.3) 

a 

In order to introduce the thermodynamic Bethe ansatz, one has to identify the string 
solutions which are supported by these equations. It is easy to see that the strings 
corresponding to (8.1) are given by 

"55 = + ^(^ + - k = l,...,n (8.4) 

""j^k — '^^r^ + ^("^ + l)/2 — i/c, k — 1, . . . ,m 

For convenience, we are going to use labels n > for it strings, n < for v strings and 
n = for the 6 solutions, which do not form strings. When going to the continuum limit, 
the equations (8.1) become equations for the density of magnon strings p„ and string 
holes Pn 

oo 

Pn + Pn — cosh TtO Sn,0 + 

(8.5) 



m=— oo 



The kernel i^n,m acts by convolution and it vanishes for nm < 0. For nm > 0, K^^^ — 
K_n-m and its Fourier transform has a relatively simple expression 

(l-Kr' = l-sI, where m = (8.6) 

and / is the incidence matrix of an A^o graph, In^m = ^n,m+i + ^n.m-i- In the rapidity 
space, the operator s can be expressed via the shift D hy s = 1/{D + D~^), where 
DfiO) = f{9 + i/2). The remaining kernels are given by 
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The thermodynamical functions can be obtained by minimizing the free energy density 
f^F/C 



/ = J del mpo cosh ne-Tj2 In (^1 + + Pn In (^1 + 



(8.8) 



Taking into account that the variation of densities of strings and holes are related by 
5pn + Spn = Z)m=-(x> ^n,mSpm, wc obtain at equilibrium 

^ cosh 7r0 Sm,o = In f 1 + ^"j + - S„,J In (l + ^) (8.9) 

^ V PraJ \ PnJ 

By denoting y„ = ^ forn 7^ and Iq = and considering the cases the m = 0, m > 0, 
m < separately, see [103] for details, one can write the above equation as 



m 
T 



cosh -kB 5^,0 = XI 

)ln(l + y„)+ln(l + y-i) . (8.10) 



Multiplying the previous equation by the operator s ^ — D + D ^, for which the l.h.s. 
is a zero mode, and taking into account that Im,n is an incidence matrix, we obtain 

In [Y^{e + t/2)Y^{e - t/2)] = In [(1 + F„+i)(l + F„_0] • (8.11) 

This is the so-called Y-system for the principal chiral field. It is interesting to note that 
the information about the potential term, which contains the size of the system, has 
disappeared from the equation and it enters to the solution only via boundary conditions. 
The equations (8.1) can be retrieved if one imposes the boundary condition at large 9 

^ g-m£cosh7r6l 5„,0 (8 12) 

The Y functions of the principal chiral field are defined on a one-dimensional lattice. 
More generally, the Y system is defined on a two dimensional lattice, Ya^g 

with some specific boundary conditions. The Y system can be reformulated as the Hirota 
equation via the transformation 

Ya,s - P^^ii^ , (8.14) 

o+ljS-' a—l,s 

where Ta^g is defined up to a gauge transformation Ta^g gi{9 + i{a + s) /2)g2{9 + i{a — 
s)/2)g^{9 - i{a + s)/2)gi{9 - i{a - s)/2)Ta,s. The equation obeyed by T is 

'^a,s'^a,s = Ta^s+lTa,s-l + Ta+l,sTa-l,s , (8.15) 
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where T^^ = Ta,s{0 ±i/2). This type of equation is obeyed by transfer matrices of spin 
chain with auxihary spaces corresponding to representations with rectangular Young 
tableaux of size (a, s) [231, 230, 232, 233]. On the border of the domain, the Hirota 
equation becomes the discrete d'Alembert equation in 1+1 dimensions, with solution 
{e.g. on the boundary a = 0) 



To,s{0)^9i{0 + is/2)g2{e-is/2) 



(8.16) 



Coming back to the Hirota system for the principal chiral field, the boundary conditions 
can be reduced, by a gauge transformation, to 



ro,.(^) = HO + is/2) , T2,s{e) = ^9 - is/2) . 
If one imposes a reality condition = Ys{u*), taking into account that 

T^^s{e + t/2)T^^,{e~i/2) 



1 + Ys{9) 



^{9 + is/2)<!?{9 -is/2) 



(8.17) 



(8.18) 



we conclude that $ and $ are complex conjugate functions^*^ , $(6*) = $*(^). The 
Hirota equation, quadratic in T, can be seen as the compatibility condition for a system 
of equations linear in T. In the simple su{2) case, this system reads, in the notation 
Tg = Tis 

Ts+i{9)Q{9 + is/2) - Ts{9 - i/2)Q{9 + i{s + 2)/2) = <l>(^ + is/2)Q{9 - i{s + 2)/2) 
Ts-i{9)Q{9 - i{s + 2)/2) - Ts{9 - i/2)Q[9 - is/2) = -^{9 - is/2)Q{9 - is/2) . 

(8.19) 

The solution of the Hirota equation can be represented in a determinant form 



Ts{9) = h{9 + is/2) 



Q{9 + i{s + l)/2) R{9 + i{s + l)/2) 
Q{9 -i{s + l)/2) R{9 -i{s + l)/2) 



.20) 



where h{9) is a periodic function, h{9 + i) = h{9) and Q, Q and R, R arc two linearly 
independent solutions of the linear system (8.19), related by the Wronskian relation 



$(^) = h{9 + i/2) 



R{9) Q{9) 
R{9 + i) Q{9 + i) 



(8.21) 



Speciahzed to the value s — 0, the equations (8.19) take the form of the Baxter equations 

Ti{9) Q{9) = To{9 i/2) Q{9 + t) + <^{9) Q{9 - t) 
Q{9 i) = Toi9 - i/2) Q[9) - '^9) Q{9) . 

The asymptotic Bethe ansatz equations of the principal chiral model (8.1) can be ob- 
tained from the Hirota equations supplemented with the boundary condition (8.12). We 
have from the definition (8.14) 



Yo{9) = 



Ti{9)T.,{9) 
<^{9)<^{9) 



g-m£cosh^0 ^ Q when £^oo. 



(8.22) 



10 



For su{n) one can take as symmetry condition T'a,s(^) = ^"n-a s(^)- 
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therefore one of T_i{0) or Ti{9) is exponentially small. Since the two objects play 
symmetric roles, one deduces that the exchange between the cases T^i{0) — > and 
Ti[6) — )■ is done by a gauge transformation which also depends exponentially in C 
Following [103], we denote the transfer matrix in the two gauges T^{0) and T^{0), such 
that 

T!i(^)->0, T^{e)^0. (8.23) 
The gauge transformation preserving the reality condition of Tk are given by 

T:(e) = g{e + zs/2yg{e - is/2) T:{e) (8.24) 
$'^(^) = g{e + i/2)g{e - i/2) (^"{6) 

^\e)^g{e + i/2)g{e-z/2) ^%e) 

so that 

T^,{0) = g{e - i/2)g{9 + i/2)Tl,{9) , (8.25) 

where g{9) is some function to be determined. Let us choose the first of the two gauges. 
Prom the Hirota equation, we deduce that 

T^{e + i/2)T^{e - i/2) = $"(^)$"(^) (8.26) 

A solution of this equation is given by Tq{9 + i/2) = $"(^). Supposing that Tq and 
are polynomials, Tq = Ylai^ " ^a) = <Pid), and = Ylji^ ~ ^j) ^^^^ 
the solution of the equation (8.19) coincides with the first line of the Bethe ansatz 
equations (8.1). The last line of (8.1) can be obtained by using the other gauge and 
Tq — Tq — (j){9) and — Ylji^ ~ ~ '^j)- We can solve the equivalent of the equation 
(8.26) by Tq{9 — i/2) — ^'"{9). this imphes that the gauge transformation g{9) obeys 

^ III 

This equation can be solved in terms of the scattering matrix S{9) = YYa^^i^ ~ ^o) 
which obeys the crossing condition 



5(» + i/W-i/2) = ^^. (8.28) 

The solution of the difference equation (8.27) is defined up to a zero mode, satisfying 
gQ{9 + i/2)gQ{9 — i/2) = 1. To fix the zero mode we use the boundary condition (8.22), 
which implies \go{9 - i/2)\'^ = g-^^^^^i^^'', and \go{9)\'^ = 1. We conclude that 

The explicit form of the gauge transformation allows to derive the central equation in 
(8.1). It suffices to consider the condition 

-1 = Yo{9^ + i/2) = '^^^~^9{0 + ^/mO - ^/2)\e=e,+./2 • (8.30) 
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and use the equations (8.19) for s — 0, the first in the u gauge and the second in the v 
gauge. 

Vacuum solution at finite C. The simplest solution of the Hirota equations for s > 
is T^__i = s and T^'^i = ^ for s < 0, which corresponds to the vacuum with strictly 
infinite C. However, T^^ should not be zero, but T^^ = TZ^\g{e - i/2)\^ ~ g-m£cosh7re_ 
Since in the vacuum T" = T^^, we are going to consider only the u gauge and drop the 
gauge index in the following. An improved solution is 

T,_i{e) = s + G{e -is/2) -G{e + ts/2) , (8.31) 

$(^) = i + G{e + io)-G{e + i) , 
$(^) = i + G{e -i) -G{e -io) , 

where G{9) is the resolvent associated with the "density" T_i{0) 

and therefore it is analytic everywhere except of a cut on the real axis. In particular, we 
have 

ro(^) ^i + G{e- i/2) - G{e + i/2) = ^{e - i/2) = l>(^ + i/2) . (8.33) 

The function Tq{9) is analytic between the cuts at 9 — ±i/2 and the functions $ and 
l» are, up to a shift, the determination of Tq beyond these cuts. The ansatz (8.32) is 
motivated by the fact that jump in Tq should be equal to TLi, c/. the first equation 
(8.19) at s = -1 and Q = 1, 

r_i(^) =To(^ + i/2) -$(^) . (8.34) 

The solution (8.31) is determined by an unique function, f{9) = T^i{9). In turn, this 
function is constrained by the gluing of the two wings in the Hirota equation, realized 
via the function Yq 

y(n. _ T,{9)T_,{9) 
Yoi9 + i/2)Yoi9-i/2) = {l + Ymr=(^P^^^§M^Y 



^9 + i/2)^9 -i/2) 
Eliminating Yq from the two above equations, one obtains formally 

($(^ + i/2)^{9 - i/2))^^^^ 

This equation can be solved numerically for arbitrary length, taking as initial vahics 
Ti{9) = 2, $(6*) = $(6*) = 1 and computing iteratively T_i{9) until it converges to its 
actual value. The results for the vacuum energy as function of mC are given in [103]. 
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Solution at finite L in the u{\) sector. In this sector, there are excitations of the type 
6a-i but no Uj or Vj. The hnear Hirota equation is again of the form (8.34), but now TQ{d) 
has zeroes on the real axis, To(6'q,) = 0. The infinite length solution is provided this time 
by Tg-i — P{6 + is/2) — P{6 — is/2), where P{0) is a polynomial satisfying 

p{9 + i/2) - p{e - 1/2) = m = Y[{d - d^) . 

a 

Let us now define the function 



with the property that on the real axis G{e - iO) - G{e + iO) = T_i(^) and G{Q) + G{Q) 
is analytic. A solution of the Hirota equation can be generated from 

To(^) = 0(^) + G{Q - i/2) + G{Q + i/2) , (8.39) 
^{Q) = 0(^ + i/2) + ^(e + iO) + ^(e + i) , 
^{B) = (^{Q-%/2) + G{Q-%) + G{B-%^) . 

The solution can be put in a determinant from via 

T^{Q) = R{Q^i{k^\)/2)-R{B -i{k^\)/2) , (8.40) 
R{d) = P{d) + ^ {G{Q + i/2) + G{Q - i/2)] . (8.41) 

The function R(Q) is a solution to the system (8.19) linearly independent from Q{Q) and 
the transition between the two solutions is done by the gauge function g{Q) 

\^Q\Q)^g{Q-i/2)R'{Q). (8.42) 

The cqTiation determining T_i recursively is the same as (8.37), and it can be considered 
as an equation for g{Q). The generic case when both Q^{0) and Q^{9) are non-trivial can 
be treated similarly, by relating the jump of the function Tq{9) at ^ ± i/2 to the gauge 
function g{9) and deriving a close equation for this function. 



8.2 Y-system for the AdS/CFT Bethe ansatz 

The Beisert-Staudacher equations (6.52) resemble to the asymptotic Bethe ansatz equa- 
tions for the principal chiral model (8.1), in the sense that the momentum carrying node 
is the central node and that the two wings are symmetric and they correspond to an 
inhomogeneous spin chain with su{2\2), respectively with su{2) symmetry. It is therefore 
tempting to conjecture that the Y-system (and therefore the T-system) corresponding 
to the Beisert-Staudacher equations is that of two su{2\2) chains, glued together via 
the central node. This hypothesis was confirmed by analyzing the bound states of the 
Beisert-Staudacher equations and writing the corresponding Y-system [246, 247]. Al- 
though the resemblance with the principal chiral field is important, there are a few 
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points which are shghtly different and which render the analysis of the equations con- 
siderably more complicated. The kernels depend on the Jukovski transform x of the 
spectral parameter, u/ g = x + 1/x. The existence of this parameter can be traced back 
to the existence of the central extension of the su{2\2) algebra [53]. The so-called mirror 
transformation [102], which transforms {p,E) {iE,ip) and which corresponds in the 
principal chiral model to the translation u — > u+i/2 is, for the AdS/CFT system realized 
by {x'^jX") — >■ {x'^, l/x~). In terms of the elliptic parametrization (6.46), this transfor- 
mation is equivalent to the translation of the elliptic argument s — )■ s -|- K — iK', which 
exchanges the real and imaginary parts of the momentum, (p, /?) — )■ (p = — f3 = ip). 
In the Jukovski plane, parametrized by the variable x, the mirror transformation sends 
the real line into the unit circle and conversely. The mirror system is not equivalent to the 
original one, since the corresponding dispersion relations are different. Let us remind that 
the particle-antiparticle transformation appearing in the crossing equation is the square 
of the mirror transformation^^ s ^ s + 2K — 2iK' and it takes {x~^, x~) — >■ (l/x"*", l/x~). 

8.2.1 Hirota equation and the Backlund transform for supersymmetric spin 
chains 

Let us first review the Hirota equation for a spin chain with GL{K\M) symmetry, fol- 
lowing [249, 250] [251]. The case relevant for AdS/CFT is X = M = 2. The domain 
where T^i^{u) ^ is non-zero is the "fat hook" with 0<a<irorO<s<Mand 
a, s non-negative. On the boundary a — the Hirota equation becomes the discrete 
d'Alembert equation, with the solution 

T^^^'iu) ^g-iu- is/2)g+iu + is/2) . (8.43) 

We can make use of the gauge freedom in order to fix the boundary condition as 

T^:^{n) = QK,M(n - ^^2) , T^6^(u) = QkM^ + ^«/2) ■ (8.44) 

For spin chains, Qk,m{u) is a polynomial. The linear auxiliary problem, which is a 
higher rank analogue of (8.19), relates the transfer matrices whose rank differ by one 
unit, (A;,m) — >■ (/c — l,m) or (A;,m) — >■ (A;,m — 1). For more details on this relation, 
called Backlund transformation, see [249, 250]. The Backlund transformation is related 
to the nested Bethe ansatz, where the diagonalization of the transfer matrix is done 
by a succession of operations which reduce the rank of the algebra by one unit up to 
rank one. For Lie super algebras, the reduction to rank one can be done in different 
ways, depending on the order in which of the two elementary Backlund transformation 
are taken. Each path connecting the point {K, M) to (0, 0) by going either to left or 
downwards corresponds to a different Dynkin diagram of the Lie superalgebra and is 
associated to a different form of the nested Bethe ansatz equations. 

The Backlund transformations can be formulated in a rather intuitive way in terms 

^^The transformation s — s + 2K has no effect on (a;+, 
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(k,m-l) 



(k,m) 



(k--l,m-l) 



(k-l,m) 



Figure 6: Two different paths for successive Backlund transformations. The blue and 
red pieces of path correspond to the two pieces of Dynkin diagrams represented in fig. 3. 



of generating functions. 



s>0 



Qk,miu) 



a>0 



TtT^u-t{a + l)/2) 

Qk,m{u - i) 



(8.45) 
(8.46) 



where again D = e*^"/^. By setting 



Yk,m{u) 



Qk,m{u + i)Qk-l,m{u - i) 
Qk,m{u)Qk-l,m{u) 

Qk,m-l{u + i)Qk,m{u - i) 
Qk,m-l{u)Qk,m{u) 



(8.47) 



one can prove the following recursion relations which encode the two kinds of Backlund 
transformation 



:i-Yk,^+,{u)D'')Wk,M 



(8.48) 



Going from the point {k,m) to {k — l,m — 1) on two different ways as in figure 6 we 
obtain 



m—l 



1^1 — -^fc,m-l-D^)(l — Yk^mD^) ^Wk,m 



k,m 



which implies the zero-curvature condition 

(1 - Yk-i,^D^)il - Xk,m-iD^) = (1 - Xk,mD^){l - Yk,^D^) 
This relation gives for the term linear in 

^m—li'^) — -^k.mi'^) ^,m('^) • 



(8.49) 
(8.50) 



^.51) 
^.52) 
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Figure 7: The T hook corresponding to the Hirota equation for psu{2,2\4). The two 
wings correspond to two spin chains with symmetry s-u(2|2), glued together on the black 
nodes. The blue (gray) path connecting the points (2, 2) and (—2, 2) to (0, 0) represent 
a chain of Backlund transformations associated to the second Dynkin diagram in fig. 1. 

0— OHg) 

Figure 8: A Dynkin diagram for a spin chain with sm(2|2) symmetry. 



which translates into relations for the Q polynomials which are equivalent to the Bethe 
ansatz equations. The zero curvature condition (8.51) insures that the Bethe ansatz 
equations corresponding to different paths (or different Dynkin diagrams) are equivalent, 
and it is related to the duality transformations described in section (2.7). From (8.45) 
we can obtain the transfer matrix for the fundamental representation 

rrK,M/ \ 

TT ~ / (Xk,m — Xk^m) (8.53) 

^^'^^^) Si: 

where Yfc^m is associated to a vertical step and ,„ to a horizontal step starting at 
{k,m). For example, for the sm(2|2) spin chain with Dynkin diagram from figure 8 , we 
have 

W2,2 = (1 - Y2,2D'){1 - X2,iD^)-\l - X,^,Dy\l - Yo^,D')Wo,o (8.54) 

so that 

T'f = g2,2(Xi,i+X2,l- 1^2,2 ->o,i) (8.55) 

^ f Qi,iQo,i I Q2,iQi,i Q2,iQ2,2 Qo,oQo,i 

— V2,2 -p^ 7^ r 



Qi,iQo,i Q2,iQi,i Q2,iQ2,2 Qo,oQo,i 

where and indices ±± represent a shift by ±i. The Bethe ansatz equations are obtained 
as the conditions of regularity for T^f{u) at the zeroes of the polynomials Qk,m{u)- 
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With the identification (5o,o = 1, Qo,i = Qi-, Qi,i = Q2, Qt,! = Qs and (52,2 = Q4, 
Qn — Y[f=ii''^ ~ '^n,j) they are equivalent with the su{2\2) wings of the one-loop system 
i.e. the first three or last three equations in (2.76). It is interesting to note that for 
K = M, if T^^^{u) is a solution of the Hirota equation with boundary condition (8.44) 
defined by Qk.k{u) = Y\j{,u — uk^k-j), then T^^^*{u) is also a solution with similar 
boundary condition defined by Qx^xi"^) — Ylji"^ ~ ''^*K,K;j)- 

8.2.2 Hirota equation and the AdS/CFT integrable system 

The results from the previous section are vahd at one loop in gauge theory. The all-loop 
dependence of the Bethe ansatz can be incorporated formally by redefining in equation 
(8.55) Q2,2{u) = R-(+\u), Ql^iu) = R-^-\u), Qo,o = B+(-\u), Q++{u) = B+(+\u) 
with 

We have Q^(u) = (-g)^-Rt\u)Bt\u) and 

Let us note that Ti i enjoys the interesting transformation property 

The Y-system for the full Beisert-Staudacher system can be obtained by gluing the T- 
functions for two su{2\2) spin chains, T' and via a gluing condition similar to (8.22), 

_ f xju- ta/2) Y f{u- ia/2) , 



x{u + ia/2)J f{u + ia/2) 
where the first two factors are zero modes, i.e. solution of the equation j-^" {" — = 1. The 

' ^ Ja + lja-l 

dependence on L is fixed by requiring a behavior of Yafl{u) at large u compatible with the 
asymptotic Bethe ansatz. The function f{u) can be determined from the requirement 
that the asymptotic Bethe Ansatz is reproduced from the equation Yi 0(^4,^) = —1. 
The energy of a state is given by 

E-Eo=j2 ^1(^4,) + E y ^ ^ ^"^^ + ^^''^""^^ ' ^^-^^^ 

where Eq is the classical dimension, = 2ig{l/x"''^ — l/x"~) and the tilde indicates that 
the corresponding quantities are evaluated in the mirror dynamics. 

There are four different presentations of the Beisert-Staudacher equations correspond- 
ing to the values 771,772 = ^1 in [52]. These different presentations can be traced back 
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to the choice T'"'' = Ta,i{x) or T^^^{x) respectively in (8.59) and they lead to different 
values for the function f^/f^. The choice r^i = 772 = 1 corresponds to a central node of 
su{2) type, r/i = 772 = —1 corresponds to s/(2) and 771772 = —1 corresponds to su{l\l). 
The asymptotic Bethe ansatz equations with sl{2) central node 



r ^ ( S{x)V Q++ Bt ^3- 

/+ \r+^+)) Q- B-,Bti BTrB 



+ ' 

3r 



(8.61) 



where S^{x) = Yif=i'^^i^j^4^,j) used the following convention for the indices of 

the nodes of the nested Bethe ansatz (1, 2, 3, 4, 5, 6, 7) — )■ (Ir, 2r, 3r, , 3/, 21, 11). 

The symmetry condition^^ ¥10(1/0;) = Y{f^{x) translates into the equation for the 
scattering factor 

S*{l/x) B+(+)R-(+) B+(+)R+(-) ' ^ ' 

which is a consequence of the crossing equation for the elementary dressing factor a{x, y) 
(6.94) 

, , , , y~ x~ — 1 1 x~ — y~ 
a{x, y)a{l x, y) = ^ ^ \' ^ _ ■ 8.63 
y+ x^ — y+ 1/X+ — y 

if we take into account that the dressing factor is a pure phase for particles (or anti- 
particles) with physical dynamics, S*{x) — S~^{x) and if we use the momentum condition 

8.2.3 Four loop Konishi operator and the Y-system 

The Konishi operator is the simplest operator which is not protected by supersymmetry, 
Tr$^. It belongs to a super- multiplet and we find it more convenient to work with a 
representative situated within the su{2) sector, Tr ZXZX. The length of the spin chain 
is L = 4 and the magnon number is M = ^4 = 2 corresponding to two roots on the 
central node, ui^4 — —U2,4 — l/2y/3 + 0{g'^). The Y functions associated to the central 
node are given by 

^«>o(x)=(^)'pn„T- (8.64) 

where the index denotes a shift of the argument by ±ia/2. The Beisert-Staudacher 
equations in the su{2) presentation are obtained if 



L-sHx)'^(?^\' SlIIlSl^ (8 65) 



where we have made a gauge transform Ti^a Ti^a/ R°'~^~^^ such that To,a = 1. Because 
we are in the su{2) sector, the B, R, Q polynomials associated to roots other than the 
central root are equal to 1. In order to compute the finite length corrections to the energy. 



^^The complex conjugation in Y*{x) is understood to send a;^ — >■ and of course it assumes g is 
real. 
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the Y functions should be computed in the mirror dynamics for the main argument, so 
we have to transform (x""*", x"'~) — >■ (x""^, For small g we have 



-9'' (8.66) 



which is in agreement with the fact that the wrapping corrections are of order g"^^ for a 
chain of length L. In particular, wc notice that for L = 4 at fourth loop order the main 
correction comes from this factor and it is sufficient to consider the other quantities at 
leading order. Using (8.45) we obtain for the transfer matrices 



rri* 

where at leading order in the mirror dynamics 



(8.67) 



(a+1)- 



with = n, =1 2i''^ ~ "^^^4^ ia/2). After simplification we have 



(8.68) 



24o(-l + 3«^ + 12(/-) 
Tia = , u ■ 8-69 



(8.70) 



The scalar factor gives, again at leading order in g 

f- _ Ql r- _ Qt'^'Qt'^^ 

Putting all the pieces together we obtain 

= / (?Mz±t^4±^) ' ^ ^ (8.71) 

where = 144 Q^'^-^)" Q("-^)+ = 16 - 48a + 60a^ - 36a^ + 9a^ + 48^^ - 1440^2 + 
72a?u'^ + 144-u^. Inserting this result in the expression of the energy (8.60) and taking 
the integrals, for example using the residue method, we obtain the wrapping correction 
(^^Konishi = 5^(648 + 864C(3) - 1440C(5)) + 0{g^^). 

At five loop order, the Bethe roots acquire wrapping corrections, which will manifest 
themselves in the individual magnon energies, that is in the first term in r.h.s. of equation 
(8.60). The second term will receive corrections from the dressing phase, since in the 
mirror dynamics, S'^{x) = l + 0{g'^). The difficult part of the computation is to solve the 
Bethe Ansatz equations o(M4,fe) = —1 beyond the asymptotic regime. This was done 
numerically in [105] and analytically in [106] and the results agree with the computation 
of [98] using Liischer techniques. 

The success of the Y system in reproducing various finite size corrections is encour- 
aging and it is likely that these equations are exact at any size. It is important to 
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understand the finite size corrections from the spin chain point of view, with a formula- 
tion of the dilatation operator as a lattice Hamiltonian which would be valid for any value 
of the coupling constant and for any length. In this direction, it would be interesting to 
clarify the connections with the Hubbard model and how to treat the length-changing 
operations for finite chains, which were elegantly circumvented by Beisert [53] in the case 
of infinite chains. 

9 Recent developments and open problems 

Another important question which awaits to be solved is that of the origin, and the 
formal proof, of integrability in the A/" = 4 gauge theory. It is commonly believed that 
the all-loop integrability of the spectral problem discussed above is due to an interplay 
between planarity and supersymmetry. The discovery [28] of another duality between 
the J\f — 6 supersymmetric Chern-Simons theory in three dimensions and strings in the 
AdS4^ X CP"' background which is presumably also integrable at all loops [30] [252] [253] 
[254], reinforced the opinion that at the core of the AdS/CFT correspondence stays the 
same tightly constrained symmetry which includes supersymmetry. 

Most of the effort concerning integrability in A/" = 4 gauge theory was concentrated 
on the problem of computing the spectrum of conformal dimensions, but there are clear 
indications that the integrability of the dilatation operator is related to other special 
properties of the theory. It was discovered [107] that the four-gluon amplitudes possess 
a recursive structure where the one-loop contribution is the main building block for 
the higher loop amplitudes. This structure led to the conjecture by Bern, Dixon and 
Smirnov (EDS) [108] that the all loop amphtude can be obtained by exponentiating the 
one-loop result. The BDS conjecture [108] was verified at strong coupling [109] for four 
gluon amplitudes. For six and more gluons, the situation is more subtle and the BDS 
conjecture does not hold anymore [110, 111]. 

The integrals entering the gluon amplitudes have unexpected conformal properties 
[112] when expressed in the variables Xi defined by pi — Xi — Xj+i where pi is the four- 
momentum of the i*'* gluon. This dual conformal symmetry, which was shown [255] to 
extend to a dual S'Uj^er-conformal theory prompted the discovery [256] of a duality be- 
tween the multi-gluon amplitudes and the Wilson loops formed by light-like segments 
with cusps, and where Xi are interpreted as the positions of the cusps. The (super) con- 
formal symmetry of these Wilson loops is broken by the cusps, and the anomalous Ward 
identities associated with this symmetry were used [257] to compute the the dependence 
of the Wilson loops on the positions of the cusps Xj. 

At strong coupling, the dual conformal symmetry and the duality with the Wilson 
loops with cusps arose via a T-duality [109], which was extended to a fermionic T- 
duahty in [113]. The action of this duality on the integrable structure of the AdS^ x 
superstring was analyzed in [114] and reviewed in [115]. 

In a recent paper [116] it was argued that the generators of the super-conformal 
symmetry and the ones of the dual super-conformal symmetry constitute the generators 
of a Yangian. This would imply the existence of an integrable structure at the level of the 
amplitudes, and not only of the conformal dimensions. Such an integrable structure was 
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identified by Lipatov [117] in a certain high energy approximation, where the amphtudes 
obey BFKL-hke [258, 259] equations. The corresponding integrable Hamiltonian is very 
similar to the one discovered in the high energy hmit of large N QCD [260, 261]. It 
is conceivable that this Hamiltonian could be generalized to a integrable Hamiltonian 
which would describe the amplitudes in the A/" = 4 gauge theory for arbitrary values of 
the momenta. 

Very recently, the problem of computing the strong coupling limit of gluon amplitudes 
with arbitrary number of legs was mapped to the problem of minimization of the free 
energy of a Hitchin system [262] . In [263] the problem was formulated in terms of a Y 
system for the Hitchin equation. 
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A Dynkin diagrams 

The commutation relations of a Lie (super) algebra q can be encoded into a Dynkin 
diagram. One starts form a commuting ensemble of generators, the Cartan generators, 
satisfying 

[Hi, Hj]^0, i, J = 1, . . . , rank(0). (A.l) 

The commutations relations of the Cartan generators with the raising/lowering operators 
can be put on the form 

[Hi, Jr] = Ti Jr (A.2) 

where r is now an index labeling the raising/lowering operator (the root index) and 
is a real coefficient (Dynkin label). The roots can be seen as vectors in a linear space of 
dimension rank(0) with basis vectors e,, r = r^ei, endowed with a scalar product. The 
non-zero roots can be divided into positive and negative roots with respect to an arbitrary 
hyperplane passing through the origin and not containing any of the roots. Among the 
positive roots, one can choose a set of simple roots as the minimum set generating the 
other positive roots by linear combinations with positive coefficients. Finally, the Cartan 
matrix encodes the scalar products of the simple roots when the square of the largest 
root is normalized to two 

Aij — Ti ■ Tj . (A. 3) 

For su{N) algebras, the simplest way to construct the root system is to introduce an 
extra Cartan element associated to the identity and extend the algebra to u{N). The 
basis of the root system can be defined using the vectors e^, i — 1, . . . , N with the scalar 
product 

ei ■ ej = Sij (A.4) 



91 



Figure 9: The distinguished Dynkin diagram for psu{2, 2|4). 
The positive roots are given by 

Si-Sj , i> j (A.5) 

while the simple roots are 

U^Ci+i-ei, i = l,...,N-l. (A.6) 

For the Lie algebras, all the possible choices of the set of simple roots are equivalent. 
The Cartan matrices and the Dynkin diagrams corresponding to the different choices of 
simple roots arc the same. This is not the case for the graded (or super) Lie algebras. 
The simplest case is that of the sl{N, M) algebras. The root space of this super algebra 
can be constructed using i — 1, . . . ,N and 5fc, k — 1, . . . , M with scalar products 

• ej = 5ij , 6k- 6i = -6ki and • 5jk = . (A.7) 

such that positive roots are given by 

Q - ej , Sk-e-j , Sk-Si , for i> j , k> I (A.8) 

The positive roots are not equivalent anymore, in particular there are roots of zero norm, 
which correspond to the odd or fermionic generators. Therefore the choice of the simple 
roots is not unique and the Cartan matrix and the Dynkin diagram are not unique. A 
particular system of simple roots is the one with the minimum number of odd roots. 
This is called the distinguished root system 

n = Q+i - Q , (^1 - Cat , Sk+i - 5k . (A.9) 

and the corresponding Dynkin diagram has only one gray (or sermonic) node which 
separates (N ? 1) + (M ? 1) white (or bosonic) node, as for example in Fig. 9. Other, 
inequivalcnt Dynkin diagrams can be obtained by permuting Sk with e^. This operation 
corresponds to a Weyl reflection with respect to an odd root. 

B Bethe ansatz equations as integral equations 

Let us consider the thermodynamical limit of the Bethe ansatz equation for a state with 
the maximum number of magnons, like the antiferromagnetic state in a su{2) spin chain 
or the fermionic state in the sii(l|l) sector. The discrete equations for the M magnons 
are 

M 

^iPkL _ JJ" ^i<f{uk,Ui) ^ (B.l) 
1=1 
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or by taking the logarithm of both sides we obtain 

M 

27r L ' 27rL 



k 1 ^ 



where we have taken as the mode numbers of magnons Uk — k. When the length of the 
chain is very large, we can take the continuum limit. We introduce the variable t = k/L 
and the density of rapidities p{u) — —dt/du. The derivative with respect to u of the 
equation (B.2) gives 

p{u) = + du' K(u, u') p{u') , (B.3) 

where the kernel is the derivative of the scattering phase 

K{u,u') = +^^^{u,u'). (B.4) 

Solving the integral equation (B.3) we obtain the magnon density for the state with the 
maximal number of magnon, let us call it the AF state in analogy with the su{2) case. 

At this point, we want to introduce vacancies in the magnon distribution, vacancies 
which will be called holes. Let us introduce a single hole at kh, such that the magnon 
mode numbers are now Uk — k + 9{k — kk), and call the new magnon density, in the 
presence of the hole, Pm{u). The new magnon density obeys 

1 1 dp 

Pm{u) + -5{u - Uh) = + j du' K{u, u') pm{u') . (B.5) 

If we denote Ph{u) = j^6{u — Uh) and parametrize 

Pm{u) + Ph{y) = p{u) + 2^h{u; Uh) (B.6) 
the perturbation Kh{u;uh) obeys the integral equation 

/oo 
du K{u,u') Kh{u';uh) . (B.7) 
■oo 

Moreover, for an arbitrary number of holes with density Ph{u) — X]n'^('^ ~ '^h,n) the 
equation (B.6) becomes 

/•oo 
du' Kh{u;u') ph{u') . (B.8) 
-oo 

which ressembles to the integral equation (B.5), except that it is an integral equation 
for the hole density. This allows to interpret Kh{u;Uh) as a hole scattering kernel (and 
to replace the semicolon by a comma). The equation (B.7) allows to express the hole 
kernel from the magnon kernel as 

Kh = . (B.9) 
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In the su(2) case the magnon kernel is a difference kernel 

1 1 

71 (u — U')'^ + 1 



K{u, u') = K{u - u') = --— — ^ , ^ or K{t) = -e"!*! (B.IO) 



where K{t) = j du e^^^K{u) is the Fourier transform of the magnon kernel. From (B.9) 
we deduce that 

e-|t| I foo g-t/2 

Kh(t) = rrr OT KJu) = — / dt COS tu . (B.ll) 

^ l + e-l*l ''^ ^ 27r Jo cosht/2 ^ ^ 

giving for the hole (spinon) scattering phase 

f°° dt e~*/^ 

The expression of the hole momentum can be also read off (B.8) 

Phiu) — 27r [ du' p(u)^=^ [ — = arctan sinh TTM + po (B.13) 

J J COShTTW' 

du' em(u') Kh(u',u) = (B.14) 

cosh TIU 

where the magnon energy for the ferromagnetic XXX model is taken to be e^iu) = 
It is interesting that by this procedure we have obtained the crossing-symmetric 
phase factor So{9) in the su{2) principal chiral model [190]. The reason is presumably 
that the low energy excitations around the antiferromagnetic state, the spinous with 
u ^ 1, arc described by a relativistic field theory with su{2) symmetry, therefore their 
scattering matrix should obey crossing symmetry. 

Finally let us specify the scattering data for the BDS chain 

K^-^it) = -e-l*l/^^ - 1 , Knt) - \ . (B.15) 

In the large g limit, one retrieves the scattering data for the Haldane-Shastry spin chain. 
If one defines the statistical parameter X — 1 — K , then one obtains as expected 

-^magnon 2 , -^spinon • (B.16) 

The generalization to higher rank cases is straightforward and we will exemplify it 
on su{2)). The equations for the two type of roots are 

pI + pI = s-rK''pl-rK'^pl (B.17) 

pI + pI - k^'pI + k^'pI 



with s the source term and the kernel matrix in the Fourier representation 
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In order to obtain the hole scattering kernel, we have to invert the equation (B.17), such 
that 

In the case of purely statistical interaction, these matrices are given by 

We deduce that the scattering kernel for the momentum-carrying holes is given by 

-2t poo -t 
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